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ABSTRACT 


An  automatic -throttle  compensation  system  has  the  capability 
of  removing  the  problem  of  aircraft  speed  instability  in  the  carrier 
landing  approach.  Various  systems  have  been  proposed  and  tested,  each 
differing  in  input  sensed  variables.  Two  representative  systems  are  in¬ 
vestigated  analytically  and  by  the  use  of  digital  programmed  Myquist 
plots  and  analog  simulations.  An  attempt  is  made  to  determine  optimum 
gain  constants  by  using  various  forcing  functions  in  the  analog  simulation. 
Comparisons  of  the  responses  of  the  two  systems  are  made  using  time 
history  analog  records.  A  digital  computer  stability  program,  applicable 
to  any  aircraft,  is  included. 
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AN  INVESTIGATION  OF  THE  EFFECT  OF  AUTO-THROTTLE  DEVICES 
ON  AIRCRAFT  CONTROL  IN  THE  CARRIER  LANDING  APPROACH 

1 .  Introduction 

The  carrier  landing  approach  speed  of  current  swept  wing  aircraft 
is  influenced  by  many  factors.  One  of  the  most  important  of  these  factors 
is  often  referred  to  as  speed  instability.  Speed  instability  can  be  ex¬ 
plained  with  the  help  of  Fig.  1 .  It  will  be  noted  that  for  Aircraft  "A"  ,  a 

Y  " 

slight  deviation  in  airspeed  from  the  normal  approach  speed  makes  very 
little  difference  in  thrust  required.  This  curve  is  typical  of  the  AF-1E 
type  aircraft.  On  the  other  hand,  a  deviation  towards  slower  speed  for 
Aircraft  "B"  causes  a  large  increase  in  thrust  required „  This  increase, 
if  uncorrected,  causes  a  loss  of  altitude.  Thus,  in  effect,  the  normal 
approach  of  Aircraft  "B”  is  made  on  the  "backside"  of  the  thrust  versus 
airspeed  curve.  This  curve  is  typical  of  the  F-8  aircraft.  Essentially 
this  steep  backside  of  the  curve  is  caused  by  a  high  value  of  C^^  , 
drag  increase  due  to  angle  of  attack  (or  lift) . 

Still  another  factor  contributing  to  airspeed  instability  is  slow 
engine  response.  With  a  comparatively  long  time  lag  in  engine  response, 
any  hesitation  on  the  part  of  the  pilot  in  correcting  an  incipient  airspeed 
loss  leads  to  wide  variations  in  altitude  or  airspeed  before  effective 
thrust  can  be  realized. 

With  an  aircraft  that  displays  speed  instability  tendencies ,  either 
a  faster  approach  speed  or  an  automatic  compensating  device  is  required 
to  prevent  undesired  altitude  loss  due  to  airspeed  deviations.  Faster 
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approach  speeds  are  not  desirable  for  carrier  approaches.  It  will  be  the 
purpose  of  this  paper  to  investigate  some  methods  of  applying  automatic 
power  compensation  to  counteract  airspeed  instability.  A  swept  wing 
fighter,,  patterned  after  the  F-8  airplane  is  used  as  the  vehicle  of  the 
investigation.  Both  digital  and  analog  computers  are  exploited  in  simu¬ 
lating  auto-throttle  and  anti-drag  system  response  to  various  inputs . 
Real-time,  pilot  controlled,  analog  simulated  mirror  approaches  are  also 
used  in  qualitatively  evaluating  each  system. 

This  investigation  is  the  first  of  a  continuing  series  to  be  conducted 
at  the  U .  S .  Naval  Postgraduate  School  to  evaluate  systems  of  automatic 
airplane  control . 
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2. 


Discussion 


The  aircraft-carrier  mirror- landing  approach  is  made  at  a  constant 
speed  and  a  constant  angle  of  glide  slope „  This  type  of  approach  allows 
precision  control  of  the  landing,  a  factor  of  the  utmost  importance  in 
carrier  operations.  However,  the  approach  is  usually  made  at  a  minimum 
airspeed  that  is  dictated  by  performance  and/or  aircraft  flying  qualities. 
Thus,  the  airspeed  must  be  monitored  closely  while  making  flight  path 
angle  adjustments.  A  simplified  equation  (given  in  Ref.  1)  shows  that 
flight  path  angle,  V  ,  is  approximated  by  the  formula: 


Eq.  (1) 


yc-  —  -  X  <^^J_ 

*  “  L  V/  JtJ 


_  3 


tJt  -  4 


4  y_/ 

W  0  ^ 


Flight  path  angle,  then,  is  dependent  on  lift/drag  (L/D)  ratio,  thrust/ 
weight  ratio,  and  rate  of  speed  change.  In  the  mirror  approach  the  pilot 
can  adjust  and  airspeed  by  manipulating  the  throttle  and  elevator.  The 
manner  in  which  each  of  these  controls  can  be  used,  however,  is  depend¬ 
ent  upon  the  L/D  ratio.  For  aircraft  with  a  high  L/D  ratio,  the  elevator  is 
used  to  control  ^  ,  while  thrust  is  used  to  maintain  constant  airspeed. 

With  a  L/D  ratio  in  the  vicinity  of  4  -  5  the  rate  of  speed  change  is  large 
during  maneuvers.  If,  in  addition,  the  L/D  ratio  decreases  with  increasing 
in  the  approach  speed  range  then  speed  instability  is  accentuated. 

A  variation  of  lift/drag  ratio  versus  lift  coefficient  for  the  F-8  air¬ 
plane  is  shown  as  Fig.  2.  It  can  be  seen  that  at  the  approach  speeds 
used,  (  1.1  -  1.2  VgL  ) »  L/D  decreases  rapidly  with  increase  in  CL„ 

As  reported  in  Ref.  1,  pilots  stated  that  in  order  to  changed  for  an 
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airplane  of  this  class,  increased  reliance  on  thrust  rather  than  elevator 
was  necessary.  In  the  case  of  the  F-8  airplane,  however,  height  response 
to  throttle  input  is  slow  due  to  small  thrust  line  angle  of  attack  and 
negligible  trim  changes  due  to  power.  The  F-8  pilot,  therefore,  is  forced 

into  using  the  stick  as  the  primary  flight  path  control.  Then,  since  the 

* 

airspeed  mode  of  oscillation  (the  phugoid)  is  lightly  damped,  the  aircraft 
will  hunt  about  a  new  equilibrium  airspeed  and  glide  slope  angle  in  res¬ 
ponse  to  an  elevator  deflection.  These  speed  changes  in  the  5  to  10 
seconds  following  elevator  input  are  disturbing  to  the  pilot. 

Another  way  of  interpreting  airspeed  instability  is  referred  to  in 
Ref.  2.  A  curve  similar  to  Fig,  1  can  be  drawn  for  thrust  required  versus 
dynamic  pressure,  q.  A  definition  of  a  stable  speed  regime  can  then  be 
postulated  as  that  portion  of  the  curve  where 


Ineq.  (2) 


To  express  this  inequality  in  terms  of  the  lift  and  drag  coefficients, 
we  use  the  fundamental  equilibrium  equations  for  level  flight: 


and 


therefore 


and 


since  W  and  S  are  constant. 
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Therefore 


UtMl 

^  d(z:) 

which  result  is  substituted  into  Ineq.  (2) 


d(t) 


>o 


d(c.) 

Since  S  is  greater  than  O,  this  implies  that 

d(Ct)/d(ci)  >° 

Taking  the  differentials  of  numerator  and  denominator: 

hi  CB-JU  *±dC,  ~  -  <V  Jcu 


rt  0K1 


—  c4  c , 


or 

Ineq.  (3)  ‘ 


oiCt> ' 


^  >  o 

Ct  dC  t 


Thus,  satisfaction  of Ineq.  (3)  implies  speed  stability. 

The  C^'  term  used  above  is  composed  of  +  C^g ,  where  C^,  is 

equivalent  propulsion  system  drag.  C.„  is  defined  as 
C  £  J.  iT  , 

Cas  -  “5  7*-  J 

or  for  constant  altitude  and  slow  speeds, 

C«  -  -  £■  iX_ 

S  S  3V* 

C^g  is  seen  to  be  proportional  to  the  negative  of  the  thrust  change 
with  respect  to  the  velocity  change.  Inequality  (3)  may  then  be  written  as 


Ineq.  (4) 


I  - 


_  Cfc-t  Cav 


ctC, 


>o 


A  more  sophisticated  approach,  based  on  the  suppression  of  altitude 
disturbances  by  the  pilot's  use  of  the  elevator,  and  originating  with  the 
basic  equations  of  motion,  would  produce  essentially  the  same  inequality, 
with  a  small  modification  factor  included.  (See  Ref.  2)  The  resulting 
inequality  is 
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£d  +  C<\S 

Cc 


However,  the  correction  term,  C  ,  ,  is  of  relatively  small  magnitude 

D/ a 

and  may  be  ignored.  From  this  same  analysis  the  time  constant, 
of  the  subsidence  of  an  airspeed  error  is  shown  to  be  proportional  to: 

Eq- (5) 

Typical  values  of  the  parameter  I  for  the  F-8  are  given  below. 

At  1 . 1  VgL,  I  =  .22  -  ,548  =  -  .328,  showing  distinct  speed  instability. 
At  1 . 2  Vg^,  I  =  o  200  -  .  222  =  -  .  022  ,  a  less  unstable  value . 

From  an  examination  of  the  speed  instability  parameter,  Ineq.  (4) , 
it  will  be  seen  that  positive  stability  can  be  produced  if  either  drag  in 

T 

trimmed  conditions  (C^q)  or  C^g  is  artificially  increased. 

CDO  could  be  increased ,  for  example ,  by  extending  the  speed 
brakes  during  an  approach.  This  would  accomplish  two  objectives;  one, 
increase  C^q  in  t'ie  stability  formula,  and  two,  cause  the  engine  to 
operate  in  a  higher  RPM  range  to  produce  the  added  thrust  necessary, 
thus  effectively  causing  a  slight  decrease  in  the  inherent  engine  time 
constant. 

Remembering  the  definition  for  C^g ,  the  drag  coefficient  corre¬ 
sponding  to  the  thrust/speed  variation,  it  will  be  seen  that  an  artificial 
means  of  increasing  thrust  with  decreasing  airspeed  accomplishes  the 
desired  objective.  This  is  exactly  what  automatic  throttle  or  approach 


d  Co 

Tel 

I  +  <~p 
9 


>  o 
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power  compensator  devices  are  intended  to  do.  Still  another  method  of 

increasing  C  is  by  the  use  of  an  anti-drag  technique.  This  may  be 
AS 

accomplished  by  having  a  drag  producing  device  (such  as  side  extending 
speed  brakes)  extended  during  the  approach.  Then,  by  using  an  airspeed 
controller  type  device  to  regulate  the  amount  of  drag  produced  (by  opening 
or  closing  the  speed  brakes)  in  response  to  thrust  (or  anti-drag)  required, 
airspeed  deviations  can  be  closely  controlled.  This  particular  method, 
within  its  limitations  of  having  enough  anti-drag  available,  theoretically 
has  the  advantage  of  a  much  shorter  time  lag  than  an  engine  response 
scheme  „ 
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3.  Artificial  Speed  Stability  through  Automatic-Throttle  Techniques 

There  is  a  choice  of  inputs  to  an  automatic-throttle  type  controller, 
including  various  combinations  of  derivatives  and  integrals  of  airspeed 
error,  pitch  attitude  error,  angle  of  attack  error,  and  normal  acceleration. 
There  are  also  many  different  philosophies  by  manufacturers  concerning 
the  ideal  combination  for  best  performance  under  all  conditions .  Two 
representative  systems ,  one  utilizing  angle  of  attack  and  normal  acceler¬ 
ation,  and  the  other  using  airspeed  error  and  pitch  angle,  were  chosen  for 
investigation . 

A)  Angle  of  Attack,  Normal  Acceleration  Inputs 

The  Specialties,  Inc.  Automatic  Power  Compensator,  (APC) ,  an 
auto-throttle  in  use  in  the  F-8  and  F-4  aircraft,  is  an  example  of  a  control 
using  these  inputs .  A  block  diagram ,  typical  of  this  type  of  system  is 
shown  in  Fig.  3.  In  this  system,  constant  angle  of  attack  is  maintained 
in  the  landing  approach  using  the  two  sensed  inputs  of  instantaneous  angle 
of  attack  and  acceleration  normal  to  the  glide  slope  in  g  units  „  The  normal 
g  input  supplies  anticipation  whenever  the  aircraft  is  not  in  a  one  g  flight 
condition.  A  thrust  command  is  calculated  in  the  system  computer  depend¬ 
ent  upon  error  signals  and  chosen  gain  and  temperature  factors.  This 
thrust  command  then  becomes  an  input  to  the  engine  which  in  turn  adds 
an  appropriate  thrust  to  the  airframe.  This  thrust,  acting  as  the  forcing 
function  to  the  airframe  modes  of  motion,  changes  the  velocity  of  the 
aircraft.  Thus  the  APC  computer  indirectly  controls  airframe  velocity 
through  the  equations 
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Eq.  (6) 


AT=  K*Ao<  +  Kn  An*  j  AV=  Kr  at 
ZW  -  K*  4K  +  K*4n* 


where  the  K  values  are  appropriate  shaping,  filtering  and  simulation 
networks . 

Or,  looking  at  it  differently,  the  C  of  Ineq.  (4)  has  now  been 

iiu 

given  a  greater  value  (  C/\ $  ~  J 


artifically,  thus  making  the  speed  stability  parameter  more  positive  or 
stable . 

The  transfer  function  equations  for  this  auto-throttle  as  supplied 
by  Specialties,  Inc. ,  are  also  shown  in  Fig.  3.  It  will  be  noted  that 
the  throttle  servo  is  simulated  by  a  .  1  second  time  constant  first  order 
lag.  Angle  of  attack  input  is  treated  by  both  an  integrator  and  lag  net¬ 
work.  Normal  g  input  is  also  integrated. 

The  results  of  the  analog  simulation  runs  using  this  type  of  control 
are  reported  in  Section  IV  of  this  paper.  Experimental  Methods  and 
Results, 

B)  Pitch  Angle ,  Airspeed  Change  Input 
Flight  path  angle ,  ,  is  given  by  the  equation 

Eq,  (7)  V  -  -&< 

During  constant  airspeed  flight,  pitch  attitude  change  (A<9")  and 
flight  path  angle  change  (Atf )  are  proportional;  but  dynamically,  A  &— 
occurs  before  A^.  Thus  the  means  of  anticipating  flight  path  angle 
changes  is  available  through  pitch  attitude  changes  recorded  by  a 
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vertical  gyro  mounted  in  the  aircraft.  This  is  the  theory  underlying  the 


pitch  angle,  airspeed  change  type  of  automatic-throttle  controller. 
Through  the  use  of  a  properly  designed  network,  these  inputs  make  it 
possible  to  minimize  transient  airspeed  changes  and  to  eliminate  steady 
state  airspeed  changes  resulting  from  variations  in  flight  path. 

The  governing  equation  for  this  controller  may  be  of  the  form: 


Here  the  tern1  is  inserted  as  a  steady  state  error  washout  term. 
The  lag  term  (  I  simulates  the  airspeed  sensor  and  filter.  The 


filter  attenuates  the  short  period  wind  gust  signals  entering  the  system 
and  acts  as  a  smoothing  influence . 

This  equation  was  reduced  to  a  linear  open  loop  transfer  function 
of  the  auto-throttle  controller.  The  block  diagram  for  this  system,  com¬ 
bined  with  the  analytically  derived  block  diagram  and  transfer  functions 
of  the  airframe  and  engine,  are  shown  in  Fig.  4.  The  stability  derivatives 
for  the  F-8  in  landing  configuration  were  used  for  the  airframe  transfer 
functions.  The  derivation  is  shown  in  Appendix  I. 

The  system  open  loop  transfer  function  was  calculated  as: 


Eq.  (9) 


LL 


U<-  |  + 
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where; 


/aT\  _  Kf^ 

i;  "  xl(i  +Ti6a)£i  TTu. jT 

/uA  _  -J££±1£IjCl±_ ^2l^=±l£l^L±  £3254}} 
UtJ/T"  (i  4  \,&uL  +  Zi.lA.x-t2-0%£+dAo24?} 

-  -  Q/4^  Cl±32J4^UlL+£l 
\StK^jA F  I&.8  (j  Ts.  ?Ul  +■  .  ~0075^F) 


 .LZ^jApr* 


AT)  -  K©- 


^7  -h  I.  I 


This  equation  was  programmed  on  the  digital  computer  to  produce  a 
series  of  Nyquist  plots,  each  utilizing  different  gain  constants. 

Fig.  5  shows  the  Nyquist  plot  for  the  equation  using  the  same  gain 
constants  as  were  used  in  the  analog  simulation  From  this  figure  the 
open  loop  gain  margin  is  seen  to  be  19.2. 


The  system  closed  loop  transfer  function,  — \  was  then  derived. 


This  is  shown  symbolically  as  follows: 


atJaf 


A  digital  program  was  devised  to  calculate  the  values  of  the  poles 
of  the  closed  loop  transfer  function  for  various  values  of  auto-throttle  K's. 
This  procedure  would  have  produced  families  of  Root  Loci  Plots.  How¬ 
ever,  because  of  lack  of  time  and  computer  malfunctions,  this  particular 
investigation  was  not  completed.  Equation  (10)  would  be  used  in  an 
optimizing  analysis  for  gain  constants . 

In  addition  the  equations  of  the  auto-throttle  wer  e  set  up  on  the 
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the  analog  computer.  This  was  incorporated  into  the  complete  system 
using  approximately  optimum,  values  of  K^.  ,  Ku,  and  Kyu.  All  simu¬ 

lation  was  done  in  real  time.  These  results  are  also  reported  in  Sec¬ 
tion  IV.  Analog  mechanization  is  shown  in  Appendix  III  „  A  composite 
block  diagram  showing  the  relative  positions  of  the  two  automatic 
throttle  systems  in  the  airspeed  loop  is  presented  as  Fig.  S. 

C)  Auto-Throttle  Controller  Linked  to  Anti-Thrust  Device 

Another  possible  approach  to  the  airspeed  instability  problem  is 
the  use  of  an  auto-throttle  controlling  an  anti-thrust  or  drag  device. 

This  could  be  accomplished  byflyingthe  approach  with  speed  brakes 
extended.  The  output  of  the  auto-throttle  would  then  reduce  or  in¬ 
crease  the  drag  by  closing  or  opening  the  speed  brakes  as  required. 
Thus,  instead  of  supplying  additional  thrust  when  the  airspeed  drops, 
the  speed  brakes  would  be  closed  a  proportionate  amount,  reducing  the 
drag  and  in  effect  accomplishing  the  same  purpose  as  an  increase  in 
thrust.  Note  that  the  time  lag  of  the  hydraulically  actuated  speed  brake 
system  is  much  less  than  that  of  the  engine,  thus  allowing  a  faster  res- 

r 

ponse.  The  time  constant  of  the  hydraulic  system  is  of  the  order  of  .1 
second,  while  that  of  the  F-8  engine  is  about  1 .16  seconds.  This 
approach  was  also  investigated  on  the  analog  computer  by  replacing  the 
engine  with  a  first  order  time  lag  of  1  seconds  (simulating  the  hydraulic 
system)  and  using  the  thrust  command  output  of  the  auto-throttle  com¬ 
puter  to  control  the  speed  brakes. 
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Obviously  the  limiting  feature  to  such  a  system  is  the  amount  of 
drag  the  speed  brakes  are  capable  of  supplying.  Of  course  it  goes 
without  saying  that  such  a  device  can  not  be  installed  on  an  aircraft 
incapable  of  completing  an  approach  while  varying  the  positioning  of 
the  speed  brakes  .  An  airplane  in  this  category  is  the  F-8  which  has  a 
lower  fuselage  speed  brake.  Insufficient  deck  clearance  is  available 
when  this  brake  is  extended,  thereby  prohibiting  its  use  during  an 
approach . 
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4. 


Experimental  Methods  ana  Results 


A)  Equations  and  Analysis 

The  basic  3-degrees  of  freedom  aircraft  longitudinal  equations  of 
motion  (in  British  notation)  used  in  the  analysis  are  shown  below.  The 
equations  are  based  on  the  small  perturbation  linear  theory.  Lower  case 
letters  represent  perturbation  quantities.  Upper  case,  subscripted  o, 
letters  represent  initial  steady  state  conditions.  Body  axes  aie  used 
throughout.  The  assumptions  made  are  as  follows: 

1.  The  airplane  has  a  longitudinal  plane  of  symmetry. 

2.  The  direction  of  the  relative  wind,  in  steady  flight,  lies  in 
the  plane  of  symmetry,  and  in  steady  state,  all  angular  velocities  are 
zero . 

1 

3.  Initial  flight  condition  is  wings  level. 

4.  Products  and  squares  of  perturbation  velocities  are  small  and 
are  ignored,  and  sines  of  all  perturbation  angles  are  approximated  by 
the  angle  itself  in  radians. 

5.  The  longitudinal  modes  of  motion  are  independent  of  the 
lateral  modes . 

6 .  Structural  deformations  are  not  considered . 

Three  degree  of  freedom  non-dimensional  equations  of  motion, 
British  notation:  (See  Ref.  3)  Figure  7  displays  angular  definitions  and 
signs . 

X  Force  Equation: 
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Z  Force  Equation; 


ur 

J  Vo 


L  4  ji/o  L  a  vr  “stju  Lkr  cfurJ  l  4 


Moment  Equation; 


where  symbols  are  defined  as  follows; 


I  = 


Wl 

fvs 


air  secs . 


dl  )_  tct(  ) 

cX.  j-  c?6t 


r 


=  -t 


£ 


u  = 

w  = 

n.  - 

AT  = 

as 

Vo 

fr  = 

u  = 
= 

VY\  = 

it  = 

c  = 
S  = 

X- 
I  = 

n= 


real  time 

air  secs 

perturbation  horizontal  velocity,  ft/sec. 
perturbation  vertical  velocity,  ft -/sec. 
elevator  deflection,  positive  down,  radians 
change  in  thrust  required  ,  lbs. 

C4  >  angle  of  attack  ,  radians 
pitch  angle ,  radians 

4  y  _ 

i  i  Inertia  Coefficient 

2 

Moment  of  inertia,  slug  -  ft 

mass ,  slugs 

tail  length,  ft. 

Y*\ 


=  relative  density 
mean  wing  chord,  ft. 
reference  Wing  area,  square  ft. 
Force  parallel  to  x  =  axis ,  lbs. 
Force  parallel  to  2  —  axis  ,  lbs  . 
Moment  about  lateral  axis  ,  lbs  . 


-  ft:. 
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h  = 
*«,  = 


XT  = 


tyn,  = 


X 


©■ 


m 


ur 


mT  = 


distance  from  thrust  line  to  CG,  ft. 


non-dimensional  Force  Stability  derivative 
and  similar  identities  obtained  by  permuting 
x ,  and  u  ,  L/s  independently . 


non-dimensional  stability  derivative  and 
similar  identities  obtained  by  permuting  x 

3- 

non -dimensional  stability  derivative  and 
similar  identify  obtained  by  permuting  x 
and  ,  T^and  T 


non-dimensional  moment  stability  derivative 
and  similar  identities  obtained  by  permuting 

u  and  Las’ 


non-dimensional  stability  derivative  and 
similar  identities  obtained  by  permuting  x, 
and 


non-dimensional  stability  derivative 


non-dimensional  stability  derivative^ h  = 
distance  from  thrust  line  to  C.G.  position 


The  necessary  stability  derivatives  were  evaluated  by  making  use 
of  a  longitudinal  dynamic  stability  Fortran  program,  designated  LONGSTAB. 
The  program  was  compiled  while  undertaking  the  course  in  dynamic  air¬ 
craft  stability  given  at  the  Postgraduate  School.  This  program  computed 
the  necessary  British  stability  derivatives  and  then  used  these  derivatives 
in  the  aerodynamic  equations  of  motion.  The  resulting  stability  quartic 
was  then  solved ,  by  the  program ,  for  the  phugoid  and  short  period 
complex  roots.  Finally,  the  periods  and  times  to  damp  to  1/2  amplitude 
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were  computed  for  both  modes.  Necessary  aircraft  parameters  used  by 
the  program  were  taken  from  F-8  data  supplied  by  Ling-Temco-Vought ,  Inc. 
or  calculated.  The  program-computed  stability-derivatives  used  in  the 
investigation  are  shown  in  Table  I,  as  are  the  other  aircraft  parameters. 

All  program-computed  derivatives,  after  necessary  conversions  were 
found  to  be  in  close  agreement  with  those  listed  in  Ref.  6^.  The  computer 
program  is  explained  in  Appendix  II .  A  sample  printout  is  shown  in 
Table  II . 

From  the  program  printout,  the  phugoid  period  is  33.5  seconds  and 
that  of  the  short  period  mode  is  6.06  seconds.  Times  to  damp  to  1/2 
amplitude  are  29.9  and  1.7  3  seconds  respectively.  Undamped  natural 
frequency,  phugoid"  *s  ca^cu^ate<^  as  radians/sec . ,  ^  as  .123. 

Short  period  Co is  calculated  as  1.19  rad. /sec.  ,  J  &s  .360. 

Ref .  6  lists  a  phugoid  6Jp  of  0 . 188  rad  /sec  ,  and  a  ^  of  .  12 . 

Short  period  results  found  were  listed  as  Cj  =1.19  rad. /sec.  ,  'S  = 

n 

0.35.  The  results  given  by  LONGSTAB  and  those  given  by  Ref.  6  are  in 
excellent  agreement. 

Thus,  as  previously  remarked,  it  is  seen  that  the  phugoid  is  very 
lightly  damped. 

*Ref .  6  reports  on  the  results  of  an  investigation  conducted  along 
similar  lines  as  the  subject  of  this  paper.  The  two  investigations  were 
performed  concurrently  and  independently  of  each  other.  After  com¬ 
pletion,  the  results  of  the  author’s  analysis  were  verified  insofar  as 
possible  with  those  given  in  Ref.  6. 
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These  non-dimensional  equations  of  motion  were  dimensionallzed 


for  the  analog  simulation.  The  engine  time  lag  and  auto-throttle  blocks 
were  next  added.  Real  time  stick  and  throttle  simulators  were  then  in¬ 
troduced  into  the  system,,  together  with  the  F-8  linearized  longitudinal 
control  system  dynamics  supplied  by  Ling-Temco-Vought .  All  analog 
runs  were  conducted  in  real  time .  The  analog  computers  used  were 
Donner  3100  and  3400  machines.  An  eight  channel  Brush  recorder  was 
used  to  record  the  time  histories. 

The  following  types  of  computer  runs  were  then  conducted: 

(a)  Transient  analysis  runs  using  step  inputs  of  elevator  deflection, 
horizontal  velocity  gusts,  and  pitch  attitude,”  first  with  basic  airframe 
alone  and  then  with  the  two  basic  types  of  auto -throttle  controllers  in¬ 
corporated.  Finally  sample  runs  were  made  with  the  speed  brake  device 
replacing  the  engine. 

(b)  Opt  imizing  runs  using  triangular  inputs  of  elevator  deflection, 
to  determine  the  optimum  gain  tor  this  type  of  forcing  function, 

(c)  Sinusoidal  gust  inputs , 

(d)  Real  time  pilot  controlled  runs  using  a  simulated  glide  slope, 
with  and  without  the  aid  of  the  auto-throttle ,  Horizontal  turbulence  was 

also  added  to  determine  the  controller  effectiveness. 

Analog  mechanization  diagrams  are  shown  in  Appendix  III .  Poten¬ 
tiometer  settings  are  shown  in  Table  1  to  Appendix  III.  Time  histories  for 
these  runs  are  displayed  in  Figs.  10  to  35. 


B)  Evaluations  of  Experimental  Analog  Runs 

The  systems  referred  to  are  as  follows; 

System  1,  Inputs  of  angle  of  attack,  and  normal  acceleration 
in  g  units . 

System  Equation:  AT  =  _J —  T  j.  _  KcvAV  ,  Ky  AM.~I 

t  +  l+.SA.  l-hy^L  J 

System  1-a  gain  constants; 

Ko<-  =  1970  lbs. /degree 

K  AV  ~  280  lbs ./knot 

K-j,  =  183  lbs. /degree  ~  sec. 

System  1-b  gain  constants; 

=  3420  lbs  ./degree 

KAv  =  72  lbs. /knot 

K-js  =  190  lbs. /degree  -  sec. 

System  2 ,  Inputs  of  pitch  angle  ,  airspeed  error 

System  Equation:  AT  —  K&-A Ktc.  ■+■ 

Gain  constants; 

Kfi-  =  318  lbs. /degree 

=  360  lbs  ./knot 

K fu,  =  22.3  lbs. /knot  -  sec. 

1)  Step  input  of  1. 32°  up  elevator  deflection 
Fig.  10  is  the  analog  time  history  response  to  a  1.3°  step  of  up 
elevator  deflection,  using  the  System  1-a  auto -throttle .  Airspeed  dis¬ 
plays  a  slight  oscillatory  tendency,  reaching  an  equilibrium  speed  oi 

+  .6%  of  V  (.8  kts  .  above  V  ).  Thrust  is  applied  smoothly  in  response 

o  ° 

to  demand . 

Fig.  11  shows  the  response  of  System  1-b.  This  system  had  the 
gain  factors  optimized  for  a  step  input  of  horizontal  gust  velocity.  Note 
the  high  airspeed  overshoot  due  to  step  elevator  input. 
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Fig.  12  is  the  same  situation  uslry?  3ys+e*rt  ?.  Airspeed  is  much 
more  highly  damped  as  compared  to  Fi  «,  Ecuil.  brium  airspeed  is 
-1%  of  Vq  (1.39  kts .  below  VQ)  „  Thrust  is  al  o  appl  ed  smoothly,. 

2)  5  knots  of  horizontal  gust  t 

Fig.  13  displays  the  record  of  -Syste  in  response  to  this 

input.  The  system  is  oscillatory  with  ar  approximate  damping  ratio  of 
$  -  .15,6^n  =  .677  rad/sec.  I^e  s'  ve  ocit  overshoot  is  61%  cf 
initial  error.  These  values  of  ca  n  co  ♦  *  ■  e  ..ess  than  optimum  lor 

this  type  of  input  disturbance. 

Fig.  14  displays  the  respo  se  or  i,a--e  ■v  „  wh  eh  has  been 
optimized  to  this  input  disturbance.  £n-ivalent  ^  is  now  „58f£^  is 
.73  rad/sec.  Thrust  response  is  sore  l  \i  oscillatory.  Final  value  of 

airspeed  is  + .  2%  of  Vq  . 

Fig.  15  is  the  record  oi  Si-  ste  i.  <,  T-  p.  n  ow  s  a  damped  response 
with  final  airspeed  of  about  -4.'  %  V  -?a*t  r  t.  »r  .es  are  less  than 

O 

that  required  by  System  1 . 

3)  4°  step  input  of  pitch  t  e 

Fig.  16  shows  the  response  o*  S ►  s*e  o  this  input.  Initial 

velocity  transient  is  +  .  8%  V  ^  or  I.)  2  ♦  .  Steady  ;  rate  error  is  + .  4%  VQ . 

Thrust  application  is  slightly  osr  .  or 

Fig.  17  shows  System  1-r  re&poi  se.  hdt'al  velocity  transient  is 

.6%  V  .  Steady  state  error  is  also  .  b  V 
o 

Fig.  18  displays  the  reopo'se  o'  3  s‘er'  .  Initial  velocity  transien 

amounts  to  +1.6%  V  .  Airspeed  <ian.ps  'o  +,4%  V4  steady  state  error  in  13 

o  ° 


seconds,  reaching  less  than  1%  error  in  9  seconds. 

Fig.  19  shows  the  response  of  the  basic  airframe  to  this  same 

input  disturbance,  without  the  aid  of  the  auto  "throttle .  Airspeed 

oscillation  is  +1%  V  ,  with  a  steady  state  value  oi  about  .5%  V  „ 

o  o 

4)  Optimizing  runs  using  triangular  elevator  input 
Fig.  20  shows  System  1-a  under  the  influence  oi  a  .02  ops 
triangular  elevator  deflection  (maximum  amplitude  1.32°)  forcing 
function.  Airspeed  error  is  kept,  within  a  range  of  +2%  to  •=■!%  of  V  . 

Thrust  varies  within  the  range  of  +1600  to  =2400  lbs.  about  the  trim 
setting.  Thrust  application  follows  elevator  movement  closely. 

Fig.  21  shows  the  reaction  of  System  I+b.  Note  how  airspeed 

fluctuates  +  4%  V  (+5.56  kts  .) .  Thrust  has  ver  y  noticeable  oscillatory 

—  o  — 

tendencies  and  gyrates  within  rather  wide  limits. 

Fig.  22  shows  System  2  under  the  influence  of  this  same  forcing 
function.  Airspeed  is  seen  to  vary  +  1.7  kts.  or  +  1.2%  of  VG.  Thrust 
varies  approximately  2400  lbs.  about  trim  thrust.  Application  of  thrust, 
however,  is  smoother  and  displays  less  hunting  characteristics  than 
shown  by  System  1-a. 

Figs.  23,  24  and  25  show  the  response  of  these  same  systems  to 
a  triangular  elevator  forcing  function  of  the  same  magnitude  as  previously , 

but  of  .06  cps  frequency. 

Fig.  26  shows  the  reaction  of  the  airplane  to  this  same  forcing 
function  without  the  benefit  of  the  auto-tlirottle  Note  airspeed  deviation 
is  less  than  that  obtained  with  System  l-b  engaged. 
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5)  Human  pilot  controlled  flight  on  simulated  glide  slope 

For  this  series  of  runs  a  pilot  flew  the  simulated  airplane  on  the 
glide  slope,  using  the  auto-throttle.,  Glide  slope  reference  was  obtained 
from  a  zero  centered  panel  meter.  The  iero  voltage  level  was  considered 
the  correct  glide  slope.  This  is  shown  on  Channel  1  of  Figs.  27,  28  and 
29 .  The  pilot  was  then  told  to  manipulate  the  stick  to  fly  the  aircraft  to 
a  new  4°  higher  glide  slope.  This  corresponded  to  a  +  14  volt  reading  on 
channel  1,  or  14  volts  on  his  reference  meter.  This  simulated  a  low 
"meatball"  on  a  mirror  approach  and  the  flight  path  correction  necessary. 
This  new  flight  path  angle  was  held  momentarily.  The  pilot  was  then  told 
to  fly  the  airplane  down  to  a  -  4°  glide  slope,  after  which  he  was  to 
return  to  the  normal  flight  path  (0  volts) . 

The  response  of  the  aircraft,  and  speed  cot  trol  is  shown  in  Figs. 
27,  28  and  29  for  Systems  1-a,  1-b  and  2  respectively. 

It  will  be  noted  that  Systems  1-a  and  2  maintain  airspeed  close  to 
within  1%  of  V  ,  whereas  System  1-b  has  much  greater  airspeed  fluctua¬ 
tions.  Thrust  application  is  again  much  more  oscillatory  in  System  1-b. 

6)  Human  pilot  controlled  flight  on  glide  slope  with  sinusoidal 
input  of  5  knots  horizontal  velocity  grusr 

Figs .  30  and  31  are  records  ot  Syste  1  runs  utilising  a  pilot 
controlled  stick  with  a  forcing  function  sinusoidal  input  of  5  knots  of 
horizontal  velocity,  frequency  .OS  cps.  The  pilot  was  again  attempting 
to  fly  the  glide  slope  4°  up  and  down. 

Fig.  32  is  the  record  of  the  response  to  this  same  forcing  function 
under  the  influence  of  the  pilot  and  the  System  2  auto-throttle . 


Fig.  33  displays  the  record  of  the  pilot  attempting  to  maintain 
original  glide  slope  while  the  aircraft  was  being  influenced  by  a  sinu  ¬ 
soidal  u  turbulence  forcing  function  ,  without  the  benefit  of  the  auto- 

throttle . 

7)  Speed  brake  anti-thrust  system 

Fig.  34  shows  the  history  of  the  speea  brake  anti-thrust  system 
in  response  to  a  5  knot  horizontal  velocity  step  input.  Mechanically, 
the  System  2  auto-throttle  was  used  to  supply  a  thrust  command,  but 
this  thrust  command  was  used  to  drive  a  1  second  time  constant  anti¬ 
drag  device,  simulating  the  movement  of  the  speed  brakes.  It  will  be 
noted  that  in  comparison  with  Fig.  15,  (standard  System  2  setup)  response 
and  damping  is  faster.  It  required  y  seconds  tor  the  transient  to  steady 

within  1%  of  V  in  the  case  of  Fig.  34,  while  it  required  11  seconds  for 
o 

System  2. 

Fig.  35  shows  the  response  of  the  anti-thrust  system  to  the  same 
forcing  function,  but  utilizing  different  gain  factors.  Response  is  more 

oscillatory. 

C)  Summary  of  Analog  Results 

From  the  foregoing  evaluations  it  can  1  e  seen  that  the  response  of 
the  different  systems  investigated  is  dependent  upon  the  sensed  variables 
and  the  applied  forcing  functions.  System  l~a,  optimized  to  an  elevator 
deflection  input,  is  marginal  to  unsatisfactory  in  performance  when 
subjected  to  a  horizontal  gust  input.  On  the  other  hand,  system  1-b, 
using  the  same  sensed  variables,  hut  with  gains  optimized  to  horizontal 
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velocity  input,  does  not  perform  well,  when  subjected  to  elevator 
deflection  inputs.  System  2,  utilizing  different  sensed  inputs,  also 
displays  varying  reactions  to  each  of  the  forcing  functions . 

In  general,  it  can  be  seen  that  System  2  has  the  capability  of  a 
smoother  application  of  thrust  in  response  to  an  airspeed  error.  System  1 
has  the  capability  of  a  more  highly  damped  response,  but  causes  a  more 
erratic  application  of  thrust.  Further  investigation  as  to  the  optimum 
values  of  gain  constants  is  required  to  accomplish  the  desired  task  more 
efficiently. 

It  is  further  recommended  that  additional  sensed  variables  or 
integrals  of  such  terms  be  added  to  the  analog  set-up  to  help  optimize 
system  response.  A  fruitful  area  might  be  to  incorporate  a  normal  g 
input  to  System  2  to  help  the  damping  characteristics.  It  is  anticipated 
that  further  investigations  along  this  line  will  be  conducted  at  the  U.  S 
Naval  Postgraduate  School.  These  investigations  will  attempt  to  utilize 
the  simplified  equations  obtained  by  the  elimination  of  the  short  period 
mode,  to  develop  a  more  tractable  analysis. 

r 

I  &  ' 
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5. 


Conclusions 


The  two  major  approach  power  compensator  systems  investigated 
appear  to  have  the  capability  of  performing  the  desired  task.  The  major 
difference  are  degree  of  anticipation  and  erraticism  of  thrust  application. 
With  proper  choice  of  gains  and  input  variables ,  it  seems  feasible  that  a 
composite  system  can  be  designed  to  provide  the  advantages  of  both 
individual  systems . 

The  speed  brake  or  anti-thrust  system  has  the  capability  of  pro 
viding  faster  response,  but  is  limited  by  drag  area  available  and  the 

physical  feasibility  of  approaching  with  speed  brakes  extended. 

The  preceding  investigation  has  firmly  convinced  the  author  that 
an  auto-throttle  mechanism  is  an  important  contribution  to  carrier  aviation 
safety.  An  aircraft  that  is  difficult  to  fly  in  the  approach  because  of 
inherent  speed  instability  problems,  can  have  its  whole  character  changed 
by  an  auto-throttle.  By  removing  the  undesired  instability  an  auto-throttle 
relieves  the  pilot  of  concentrating  on  airspeed  control,  thus  allowing  more 
precision  in  glide  slope  and  line-up  control.  This  ad  vantage  by  itself, 
should  tend  to  reduce  the  frequency  of  ramp  accidents. 

In  addition,  when  an  automatic  landing  system  is  considered,  the 
pilot-controlled  airspeed  loop  must  be  replaced  by  an  automatic  system 
to  provide  consistency  and  predictability  of  response. 

Still  another  advantage  of  the  automatic  airspeed  control  system  is 
the  reduction  of  minimum  approach  speed  l  clew  that  recommended  by  the 
pilots  when  airspeed  is  controlled  by  manual  throttle  manipulation.  In 


addition  the  automatic  device  can  he  design  e.3  to  provide 
and  instant  response  that  can  only  he  obtained  front  m  .oh 
experience.  Thus  it  is  seen  that  the  auto-t:  nettle  car*  ac 
equalizer  in  carrier  approach  trainin'/. 


the  anticipation 
actual  pilot 
h  as  *he  gre.it 


TABLE  I 


Computed  F8  Stability  Derivatives  anc:  Associated  Aircraft  Parameters 
British  Non-  Dimensional  Notation 


Xu  = 

-.2630 

Xw" 

.16996 

Xq  = 

0.0 

xe-  = 

-.450 

^u- 

-.90 

^  w  ~ 

-1.5915 

2q  ~ 

-.5878 

II 

-.0640 

M  = 
u 

.006052 

M  - 
w 

-.15896 

ESI 

II 

-.1422 

Mq  = 

-.5119 

/6i  “ 

54.466 

1P  = 

.7081 

t  = 

3.2772 

ii 

h-3 

O 

.90 

II 

Q 

O 

.263 

CT  “ 

.  195 

Sw  = 

375  ft2 

- 

-.380 

6b 

8.1° 

oVlot 

II 

.4772 
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C  position  =  24% 

Weight  =  22,000  lbs 
u  =  139  kts  =  234  ft/sec 
£>t=  -.1949 

-.349 

Xj  =  . 2046X10-4 

Mt  =  -.635X10-6 
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APPENDIX  I 


Derivation  of  Transfer  Functions 


The  dimensional  derivatives  shown  in  Table  I  were  used  to  calculate 
the  analytic  transfer  functions „  Calculated  values  of  these  derivatives 
using  F8  data,  are  shown  in  Table  II. 

The  longitudinal  Equations  of  Motion  written  in  Matrix  Laplace  Form 
are  (neglecting  negligible  terms); 

A-X*  -Xv  +3  \/W*)\  (o  Xar\ 


-2*.  4-2. 


-4.V0 


|  OjU) 

M«.  )J 


ift  0 

flse  At. 


7 


The  Transfer  function  /  \  Airframe  then  can  be  written  as 

U<d 

(in  determinant  form); 

0  "Xvr  4  3 

A-Xur 


Lc 


D, 


Where  D^  =  determinant  of  homogeneous  equations 

Expanded  Dj  =  /\aa  -+$4  +  C^L  -+D  A  +  E 
Where  A  =  1 

B  =  -  (M%.  +Xu_  +Zur  +UoM 

C  =  Mf.2u.-Ui  H„*Xu(Mf>2*r+UMir)-X*-Zu 

d  - 

e  =  ^  (  Mm-2w  -  Mu,  2^) 


*T7 


Expanded  /\lu_  —  +■  Cu.  -d.  "t  D  LA- 

Where  Bu  =  ur 

cu=  +t1i'(U*X^-3) 

Du  =  0  C Mfe "?ur  ~rlu^?rfe) 

Substituting  values  for  these  derivatives 

U\  _  ,n? i [ y <  4^-q 4  f85.il _  _ 

^Wwne  a**  .gtS*3  14*  •*•■  ■*•. 


By  a  similar  procedure  , 

/©_\  _  -i.SoC^^-.ZISS  4  +  .OOS&&) 

Oi 

(Positive  elevator  deflection  is  down) 

The  determinant  for  calculating  /-—) 

Ul/y 

XlT  ~)Cur 
O  -4  -  las  —  -4i1 l^o 
Mar  -^Alur^Mur)  <4.*- My. 


'Airframe- 


4- 


can  .  -e  written  as 


D, 


Expanded: 


6 


^U-A  _  /lr-4.^  +  Cr4  'hDr 

ATAiVfwme  h. 


Where 


At  -  X 


or 


I  (a/- 


Bx=  -X«t[  Z^+M^J  +  LTo!vli 

cT  =  X4r[i^-U>Mw]i-M6T[l/<,Xw-9] 

Dt  =  AW  3  £ur 

After  substituting  in  values:  *  , 

/U\  ^  *0ol45  +■  >8o^^lA-  "h  ~t~ \Q4-$Oj 

UuAl^Mwt  ~  £>, 
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APPENDIX  -X 

TABLE  I-A 

LONGITUDINAL  DIMENSIONAL  STABILITY  DERIVATIVE 


PARAMETERS 
{STABILITY  AXIS  SYSTEM) 


IN  TERMS  OF  BASIC  STABILITY  DERIVATIVES 

IN  TERMS  OF  NON- 

QUANTITY 

DIMENSIONAL 

NON -DIMENSIONAL 

DIMENSIONAL 
STABILITY  DERIVATIVE 

DEFINITIONS 

UNIT 

PARAMETERS 

x. 

j  i  ax 

’*  5u 

JL 

sec 

(-CD-C0a) 

■ix. 

r 

*. 

.  1  ?X 

m  cvr 
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(*'«•«) 
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.1 
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ft 

see 2  rad 

.Si  x, 

t  z» 

.1  2? 

ii 
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.1  z. 
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9  vW 
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■3P  (-%  Co) 
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.1 
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.  1  ?z 

e  Zq 

ft 
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-r  2« 

-1  ZL 

a 

I 

sec  2  rad 

.  /  .  c  \ 

^  >«.) 

■  S! 

r  ** 

> 

w. 

.x  p 

I,  iU 

1 

sec- ft 

-•£  “• 

M. 

r.  JL 

If  ^ 

1 

sec-ft 

*  M7C*a 

Bfc 

■  i  ^ 

1  f 
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ft 

c*& 

■*j  m% 

.  i  ^ 

*  I,  E<J 
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sec  rad 

.^c 

«,  % 

.In 

T  * 

1 

sec2rad 

■  e?p> c. 

2T7  •«* 

U  m 

■  —  in  r 
rc  *k 

f 
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APPENDIX  -X 

TABLE  J-3 

LONGITUDINAL  NON-DIMENSIONAL  STABILITY  DERIVATIVES 
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APPENDIX  II 


DYNAMIC  STABILITY  FORTRAN  PROGRAM 

Longitudinal  -  Program  LONGSTAB 

This  program,  once  entered  with  required  airplane  parameters, 
computes  the  longitudinal  stick  fixed  and/or  stick  free  (when  pertinent) 
stability  derivatives,  and  the  corresponding  stability  equation  coeffi¬ 
cients.  It  then  solves  these  equations  for  roots,  periods,  and  times  to 
damp  to  1/2  amplitude  for  both  the  phugoid  and  short  period  modes.  In 
addition  the  value  of  Routh's  discriminant  is  determined.  The  sequence 
of  computations  can  be  repeated  for  as  many  different  flight  regimes  as 
desired . 

There  are  a  total  of  five  data  cards  required  for  the  stick-fixed 
only  solutions ,  with  an  additional  two  cards  required  for  a  stick-free 
solution.  Thus  for  a  combination  stick-fixed  and  free  solution,  there 
will  be  a  total  of  seven  data  cards  required  for  each  flight  regime 
desired . 

The  initial  data  card  controls  the  types  of  solutions  desired 
(stick-fixed,  free,  or  both),  and  the  number  of  consecutive  runs  to  be 
made.  A  list  of  program  symbols  with  their  meanings  appear  in  Table  I 
to  this  Appendix . 

The  second  card  contains  an  alpha-nuireric  run  identification  or 
an  arbitrary  title . 

The  third  card  is  the  first  of  the  general  input  data  cards .  This 
card  should  have  entered  on  it: 
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St(tail  area) ,  Sw(wing  area)  /  cLck  t'  tail"  SB  ^body  area)  *  *B 
(body  length) ,  ^  ,  g (acceleration) . 

The  parameters  are  entered  eight  to  a  card,  each  being  allowed  10 
columns.  There  are  80  columns  on  each  card.  Each  parameter  may 
be  entered  anywhere  in  its  assigned  ten  columns,  but  must  have  a 
decimal  point  included  in  it.  This  same  format  is  followed  for  each 
of  the  succeeding  cards  „ 

The  fourth  card  should  have  the  following  aircraft  parameters 
entered  on  it: 

W  (aircraft  weight,  lbs.);  ;  0o  in  degrees;  h,  (distance  from 

2 

c.g.  to  thrust  line,  ft„);  C  (coefficient  of  thrust);  C  ;  Iw  (slug  ft  .); 

t  mty  11 

V  (forward  velocity  ft. /sec.) . 

The  fifth  card  should  have  the  following: 

CD;  CLw^;  el TA  ;  CL  wing;  C  (wing  chord);  lt  (c.g.  to  a.c.  tail,  ft.); 

Xc . g .  (x  distance  pf  c .g . ,  ft .);  XAC  (x  distance  of  wing  a.c.  in  ft .) . 

Cards  2  through  5  are  repeated  for  each  different  run ,  stick  fixed . 

If  it  is  desired  to  compute  a  stick  free  solution  simultaneously,  then 

two  additional  data  cards  (6  and  7)  are  required  for  each  run. 

Card  6  contains  the  following  parameters: 

2 

C.p  (tail  chord,  ft.);  K  (radius  of  gyration,  squared)^^;  B  ;  B  ;  B’; 

®  J,  2  2 

z>l  ;  • 

Card  7  contains  the  following  parameters: 

2 

m\  '  m\.  '  me  (mass  elevator,  slugs);  S  (area  elevator,  ft.). 
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Thus  for  one  aircraft  flight  regime  (say  Sea  Level,  Mach  0.6), 
for  both  stick  fixed  and  free ,  there  would  be  required  data  cards 
1  through  7.  For  each  succeeding  run,  cards  2  through  7  would  have 
to  be  repeated  with  the  new  parameters  entered.  Card  one  would  have 
to  be  suitably  prepared  to  reflect  the  total  number  of  runs  and  mode  of 
runs  desired.  See  Fig.  1  this  Appendix  for  sample  data  cards.  Note 
the  first  card  is  of  format  13  for  the  number  of  runs  ,  and  II  for  mode  of 
operation.  All  other  data  cards  are  of  8F10.0  format  (each  card  holds 
8  fields  of  10  columns,  each,  one  data  number  to  be  entered  per  field, 

t' 

in  floating  point  format;  i.e. ,  with  a  decimal  point  somewhere  in  the 
number.) 

The  program  prints  out  the  calculated  stability  derivatives  as 
well  as  the  computed  solutions  to  equations. 


!* 
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TABLE  I 


PROGRAM  SYMBOLS  AND  MEANINGS 

Program  LONGS  TAB 

Program 

mnenonic 

Meaning 

NOS 

Number  of  runs  to  be  made.  13  format,  card  1 

MODE 

Decision  stick  fixed  (1),  or  stick  free  (2).  11.  card  1 

ST 

St,  tail  area,  card  3. 

SW 

S.  .  wing  area,  card  3. 

CLALFAT 

,  Tail  lift  curve  slope,  card  3. 

ETATAIL 

,  Tail  efficiency,  C^/C^.  card  3. 

SB 

Sg,  Body  area  card  3. 

BL 

lg  t  Body  length  card  3 

RHO 

^  ,  Atmospheric  density,  card  3. 

G 

g,  Acceleration  due  to  gravity,  card  3. 

W 

w,  aircraft  weight,  card  4. 

EPSALFA 

card  4. 

THETA  1 

,  initial  angle  of  pitch,  degrees,  card  4. 

YH 

h,  distance  from  c.g.  to  thrust  line,  card  4. 

CT 

coefficient  of  thrust,  card  4. 

CM  ALFA 

C  ,  per  radian .  card  4 . 

mt/. 

YI 

I  ,  moment  of  inertia  .  card  4 . 

V 

V,  forward  velocity,  ft. /sec.  card  4. 

CSUBD 

Cg,  coefficient  of  drag,  card  5. 
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CLALFA 

,  winglift  curve  slope,  card  5. 

EPIA 

SjfA,  card  5  . 

CSUBL 

CLwing'  card5' 

C 

c,  wing  chord,  card  5. 

TL 

lt,  c.g.  to  a.c.  tail,  card  5. 

XCG 

x  _  ,  x  distance  of  c.g.  card  5. 
c.g. 

XAC 

xa  c  >  x  distance  of  wing  a.c.  card  5, 

CTAIL 

Cj.,  Tail  chord,  card  6. 

XKE 

2 

k  ,  radius  of  gyration  squared,  card  6. 

UE 

i  relative  elevator  density,  card  6 

B 1 

B,  ,  C,  .  card  S . 

1 

B2 

B2'Chie  •  card6- 

B2PRIME 

B„  ,  1/2  C.  .  .  card  6 . 

2  h£e 

ZETA 

/  card  6. 

ZETADOT 

card  6 . 

AMETA 

,  card  7  . 

AMETADT 

V)  *  card  7 . 

XME 

m  f  mass  elevator,  card  7. 
e 

SE 

Se,  area  elevator,  card  7. 
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F outran  Program  L0NG5T/AB 


JOE  EELL»ICNGST/E  US  HMTES 

PROGRAM  LONGS  TAB 

PROGRAM  FCR  THE  COMPUTATION  OF  THE  RCC.TS  CF  THE  CH AR ACTER I SH I C 
DYNAMIC  LCNGITUCINAL  STABILITY  (STICK  FIXEC)  ECUATICNS 
Cl  MENS  ION  KAPPA  (1  0  )  > C PR (  129  ),  CPI  M29)  , ROOTR  (  1  28  )  ,RCCT  I  (  128) 

01  MENS  ICN  CRR(129),CR  I(  125),XR(4),XI<4>  ,FXP  |4)',FXI(4)  ,SR(3)  I  St  <3  ) 
READ  t 6,  NOS,  MODE 
66  FORMAT  <13, I)  ) 

CO  777  J  *  1  *  NCS 
IC0UN1  *  C 

OCALL  INPLT1  ( C  LAI  F  A  ,  E  PI  A  ,C  »  TL  ,  XCG  ,  XAC  ,ST  *  SW.CLAU  fATt.ETiAT  A  II)  $8  »EL  « 
1  RH0,C-,W,EPSALFA,THETA1,YH,CT,CMALFA,YI,V.CSl£ECVCSUeL, KAPPA  ) 

COMPUTATION  CF  STABILITY  DERIVATIVES 
Q  *  .5  *RHO  *  V**2 
TCARAT  *  W  /( G*RHC«V*SW) 

THETAC  *  THETA  1  /  57.29577*35 
XU  *  -CSUBO 

XW  *  .5«CSU8L«( 1. C-(2.0*CLALFA/EPIA)) 

ZW  =  -.5MCLALFA  +  CSUEC  ) 


*  (  (XAC- XCG  I/O 


) 


APRIME  *  .414  +  2.0*<  (XAC-XCG)ZC) 

EPRIME  =  -.58  -  CLALFA  «  APRIME 
WINGMC  *  0.25  *  <<  C/TL  )**2  MBPRIME 
TAILMC  *  -.50*  CLALFAT  *  (ST/SW)  *  ETATAIL 
BODYMC  -  -0.C1*  <SE/SW)*((EL/TL)**2 

TOTALMQ  *  WINGMC  ♦  TAILMC  +  BOCYMQ 
Y I B  -  YI/I ( W/G  )*TL»*2 ) 

WINGZC  *  -.25»CLALFA*(C/TL)*(.44+(XAC-XCG)«2.C/C) 

TAILZC=  -.5*CLALFAT*(ST/SW)*ETATAIL 
U1  *  W/(G*RHO*SW*TL) 

EARMWCT  =  -.5*CLALFAT*ETATAIL*I ST/SW) *EPSALFA 

TOTALZQ  =  WINGZC  +  TAILZQ 

ZU  =  -CSUBL 

XTHETA  *  -.5 *CSUB  L 

ZTHETA  =  -.5*CSUBL*TANF( THETAC) 

AMU  =  -CT  *YH/TL 

AMW  =  .5  *  CM ALFA  *  C/TL 

PRINT  5C5,ST,SW,  CL  AL  F  AT  ,  ET  AT  A  I L,  SB  ,BL  ♦  RHO  ,G  ♦  W,'EPS  ALFA ,  THET  A  1,  YH  , 

1  CT,CMALFA,YI 

5050F0RMAT  ( 1 H 1  /40H  INPUT  PARAMETERS  /✓ 

1  1 1 H  TAIL  ARE A=  F8.4.2X11H  WING  AREA*  F9.4,2Xt14H  CL  ALFA  TAIL* 

2  F9.6,2X,  1CH  ETA  TAIL*  F9.6,2X,11H  BCD Y  AREA*  F9.2  / 

3  13H  ECCY  LENGTH*  FE.5.2X,  SH  RHO=  F1C.7.2X,  5H  GEE*  FE.S,2X, 

4  eH  WEIGHT*  F12  .4,  IX, 3HLBS  2X,18H  D  EPSILCN/O  ALFA*  F9.6  / 

5  2 1 H  INITIAL  PITCH  ANGLE*  F9 .6 ,  IX ,7HGEGREES  2X.22H  THRUST  LINE  Cl 
6 STANCE*  F8.5 ,  1X,4HFEET  2X,16H  THRUST  CCEF  CT*  FE.6  / 

7  eH  CMALFA*  F9.6.2X.1CH  I  SUB  YY=  F  12 . 3 , IX ♦ 1 SHSLUG  FT 

WRITE  CLTPUT  TAPE  5,  101 

1C1  FORMAT  I///  47H  STABILITY  DERIVATIVES 

WRITE  OUTPUT  TAPE  5,  1C5  ,V, TCARAT. C 

105  FORMAT  (19H  FCRWARC  VELOCITY*  F12.3,10H  T  CARAT*  F1C.6,  4H.  C* 

1  F  14.4  /  ) 

WRITE  CLTPUT  TAPE  5,  106,  CSUBL,  CSUBC 

106  FORMAT  (  9H  C  SUBL*  FlC.6,  9H  C  SLBC*  F1C.6  /) 

WRITE  OUTPUT  TAPE  5 , 1 C2 , XU, XW , XQ, XTHE TA , ZU, ZW ,TCTALZC 

1 020 FORMAT  (5H  XU*  F7.4,6H  XW*  F7.5,6H  XC=  F7.4,10H  XTHETA*  F7.5, 

1  6H  ZU*  F7.5 ,  6H  ZW*  F7.4,12K  TOTAL  ZC=  F7.4  /) 

WRITE  OUTPUT  TAPE  5 , 1 C3 , ZTH ET A , AMU, AM W , B ARM WCT,  TCT ALM C  ,U  1  , Y  I  E 
1 030 FORMA T  (  9H  ZTHETA*  F7.4,  6H  MU*  F8.6,  6H  MW*  F7.5,14H  MeAR  W  C 
1 0T*  F7.4,  6H  MQ*  F7.4,  9H  MUCNE*  F7.3.13H  I  SUB  BETA*  F7 .4  //) 

OCALL  PRERCOT  ( XU, X W, XQ, XTHETA ,ZU ,ZW, TCT AL ZCt ZTHETA . AMU, AMW, B ARM WCT 
1 ,T0TALMC»U1 ,YIB,TCARAT,A1 ,B1 »  Cl , D 1 , E 1 ,  CRItCRR  ) 

ROUTH  =  CRR(2)«  C RR 1 3 ) »  CRRI4J-  CRR(1)«  CRR (4 )«*2-CRR ( 2 )**2*CRR( 5) 
WRITE  OUTPUT  TAPE  S,  108,  ROUTH 
FORMAT  (////35H  THE  VALUE  OF  ROUTHS 


SCUAREC//) 

//) 


108 

1  1 

1C9 

12 

110 

13 


THE 


IF  IRCUTH )  11,11,12 
WRITE  OUTPUT  TAPE  5, 
FORMAT  (28H 
GO  TO  13 
WRITE  OUTPUT 
FORMAT  (26H 
CONTINUE 
‘  N  *  4 
CALL  FOLYMUL 


DESCR IMINANT 


F  10.2 


) 


1C9 
AIRCRAFT 


TAPE  5,  11C 
THE  AIRCRAFT 


IS  UNSTABLE 


IS  STABLE 


/?/  ) 


77/  ) 


ICRR  »CRI ,N, KAPPA, TCARAT) 


APPENDIX  II,  Fig.  2 
80 


) 


DECISION  FCR  LONG  CR  LATERAL,  STICK  FIXED  (999  )',OR  STICK  FREEI2C1) 
GO  TO  (777,201 ,202, 2C4) ,MCDE 
2  Cl  CONTINUE 

WRITE  CL T PUT  TAPE  5,301 

3010FCRMAT ( ////  80H  STICK  FREE  LONGITUCINAL  DYNAMIC  STABILITY  CCEFF 
1ICIEN1S  ////  ) 

WRITE  OLTPUT  TAPE  5f,  302,  KAPPA 
3C2  FORMAT  (  1 CA8  ) 

OCALL  IN  PL  T  2  (CTAIL'»XKE»-UE»B1»B2»B2PRIME»ZETA»ZETACCT,AMETA» 

1  AMETADT  ) 

XKE  =  SCRTF(XKE) 

ZQ  =  TCTALZC 
AMQ  ■  TCTALMQ 

OCALL  PRER0T2  ( CTA  ILf, EPS AL FA ,  XK E ,  Y 18 ,U  1 , UE , TL  ,  E  V, 82 ,82PRIME,ZK, 

1  ZTFETA,ZQ,ZETA,ZETADCT, AMW,AMQ,BARMWCT  ,AME<TA  ,'AMETACT  ,  CRR,CR-I 
N  =  5 

CALL  PCLYMUL  (CRR , CRI ,N, KAPPA , TCARAT ) 

777  CCNTINLE 
END  FILE  5 
ST0P2 

202  ST0P4 

203  ST0F5 

204  S  TOPE 
END 

OSUBROLTINE  INPUT1  ( CLALFA , EP I A ,C  »T L, XCG , X AC , ST , SK ,C L AL F AT , ET AT  A  IL  , 

1  SB„BL,RFO,G,W,EPSALFA,TFETAl , YH , CT ,CM AL F A , Y I , V , CSU8C ,C SUE U K A Ff A ) 
Cl  PENSION  KAPPA(IC) ,CPR(129) ,CPI (  129) ,ROCTRI128) ,RCETI(128) 
DIMENSION  CRR(  129  )  ,CR  I(  129)  ,XR (4 ) ,XI ( 4  )  ,FXP  (4  )',FXI  (4  )  »  SR  (  3  ).»  S  I!(3  ) 
READ  6  5,  KAPPA 
1C4  FORMAT  (  8F10.C  ) 

READ  10*4  ,ST,SW,CL  ALFA  T,  ET  AT  A  I  L .  SB,BL,RFO,G 
READ  1 C  *4  ,  W,EPSALFA,THETAl,YH,CT,CMALFA,YI,V 
READ  10*4, CSLBD,  CL  ALFA, EPIA,CSUBL,C,TL, XCG, XAC 
65  FORMAT  (  1CA8) 

RETURN 

END 

SUBRCLTINE  PREROCT  (XU,XW, XC , XTH ETA, ZL , Z N ,T0 T AL ZQ , Z TFE T A , AMU, AMU, 
1EARMWCT,T0TALMQ,U1,YIB, TCARAT, A1.81,C1, 01, El,  CRI,  CRR) 

DIMENSION  KAPPA<10).CPR( 129),CPI< 129 ) , RCCTR (  1  28 ) ,RCCTI(  128) 
DIMENSION  CRR(129 ),CRI ( 129) ,XR(4  ),XI( 4) ,FXR(4  ),FXI ( 4) ,  SR  (3),  SI  (3) 

A 1  =  l.C 

ZCMU  =  (l.C  + ( TCT ALZC  /  Ull) 

2  El  =  -(XU  +  ZW)-( TCTALMQ/YIB)-(ZQMU  *  EARMWCT/YIB  ) 

3CC1  *  (XL*  ZW  -  XW*  ZU  )  +  (TCTALMC/YI  E)  »  (XU*  ZW)  ♦(  XU*ZCMU  -  ZU  * 
KXQ/Ul)  -  ZTHET  A)  *  (  8  ARMWDT/Y  I  B  )  - ( U 1  * AMW/ Y 1 8  ) * < ZQMU  +  XC/U1) 
v  40D 1  =  -(TCTALMC/YIE)*  (XU*ZW  -  XW  *  ZU )  +  (-XTFETA*  ZU  ♦  Z THE T A  * 

1  XU  )  «  (  E  ARMWDT  /  Y  I E  )  ♦  ( U 1  *  AMW/ Y  IB  )  *  (  XU*  ZCMU  -ZU*  XC/U1  -ZsTHETA)* 

2  (  U  1  *  A  M  L  /  Y  I  8  )  *  ( — X  T  H  E  T  A  -XW*  ZCMU  +  ZW  *XQ/Ul  ) 

50E 1  =  (Ul*AMW/YI8) * (— XTHET A*ZU  ♦  ZTHETA*XU )-( U 1 *AKU/ Y 1 8 ) *J-XTH ET A* 
1ZW  +  ZTFETA*  XW  ) 

WRITE  CUTPUT  TAPE  5.1C1,  A 1 , 8 1 , C 1 , D 1 , El 
1 C10FCRMAT  (/3H  A*  1X,F13.5,4K  B*  1 X , F 1 3  *  5  *  4H  O  1X.F13.5,  4H  C« 

1  1X,F13.5,4H  E*  1X»F13.5  // ) 


1 


CRR ( 1  ) 

S 

A  1 

CRR ( 2  ) 

s 

B 1 

CRR ( 3  ) 

s 

Cl 

CRR ( 4  ) 

s 

D  1 

CRR (5  ) 

s 

El 

CRI (  1  ) 

s 

0.0 

CRI (2  ) 

s 

C.C 

CRI  (3  ) 

s 

C.C 

CRI (4  ) 

s 

o.c 

CRI  (5  ) 

= 

o.c 

RETURN 
END 

OSUBROUT INE  INPUT2  (CTA  IL»XKE,UE,B1  ,B2, 82  PRIME,  ZET  A , Z8T  ACCT  ,  AM.  ETA', 

1  AMETADT  ) 

SUBROUTINE  BRINGS  IN  STICK  FREE  FACTORS 
OREAC  ItCTAILt  XKE,  UE,  8  1 ,82  ♦  82 PR  IME  ,ZET  A,ZE7ACCTC,iAME  TA  ,  AMETADT  , XME  , 
1  S  E 

1.  FORMAT  (  8F10.C  /4F1C.0) 


! 


81 


WRITE  OUTPUT  TAPE 
2CF0RMAT(////1CX70H 
1 

END 


5!,  2 


INPUT 

) 


PARAMETERS  STICK  FREE 


SUBROUTINE  PREROT  2  <  C  T  AI  L  ,  EPS  AL  F  A  ,  XKE  ,  V  I B  ,UT  ,UE't  TL  ,  E  1  ,  E2 ,  E2P  P  I  ME ,  Z 
1  W, ZTHETA, ZQ, ZETA,  ZET ACOT , AMW , AM C» BARM WOT , AMET A, AME TAD T , A , E  ) 

;  SUBROUTINE  CAUCULATES  STABILITY  STICK  FREE  GUUNT  IC  COEFFICIENTS 

CIMENSICN  A ( 12  9),  B1129) 

UF ACT CR  =  2.0*  UE  /(Ul**2) 

UFACT  =  UFACTOR 
FACTOPK  =  l  XK  E  /  TL)**2 
XKFACT  ■  FACTORK 
CT FAC T  =  CTAIU  /  TU 
ZF  ACT  =  (  1  .C  +! ZC/U1  ) ) 

EPSFACT  =  1.0-  EPSAUFA 
WRITE  OUTPUT  TAPE  5,2CC1, 

2001 OFORMAT  ( 49H 

1  9H  UF  AC  TOR=  FI  4.7,  9H 


2 1  OH  ZCFACT CR=  FI H . 7 .  16H  EPSAUFA 

A  COEFFICIENT  CALCULATION 


UFACT,  XKFACT  ,CT  FACT  ,ZFAC.T*EPSF  ACT 
,  FACTOR  CALCUL AT  IONS 

KFAC TOR=  F  1U  . 7  ,  1CH  CTFACTCR*  F14.7 


// 


FACTOR*  F  1  *4 . 7 


.// 


A( 1 )  =  UF AC  T  *XKF  AC  T  *  YI8 
B  COEFFICIENT  CAUCUUATION 

OA l 2 )  =  UF  AC  TOR  *  XKFACT*!  BARMWDT  *( ZE TACCT *CTFACT/U1 -ZFACT ) -AM C-ZW * 

1  YIB  +  AMETADT*  CTFACT  )  -  B2PRIME  *  CTFACT  *  YIE  /  Ul 
C  COEFFICIENT  CAUCUUATICN 

CA(3)=UFACT  *  X  K  F  AC  T  *(  ZW  *AMC  +  AMWM  ZETACOT  *CTFAC  T-  U1*(ZFACT) 

1  )  +  EARMWDT  * ( ZE  TA  -  ZTHETA1+  Ul  *  AMETA  -  AMETADT  *  CTFACT  *ZW  ) 

2  +  E2FRIME  *  CTFACT  /  Ul*<  AMC  +  ZW*YIE  +  ZFACT  *  BARMWDT  )  -  E2  * 

3  YIB  -B 1  *  l  AMETADT  *CTFACT/Ul  ♦  ZETADCT  *CTFACT/0'1  »l  YIE*EPSFACT 

4  +  BARMWCT  /Ul)) 

D  COEFFICIENT  CAUCUUATION 

CA14)  =  UFACTOR  *  XKFACT  *(  U 1  * AMW « ( ZE T A-ZTH E T A  )-  U1*AMETA«ZW  ) 

1  -B2PP I  ME  *  CTFACT  /Ul  *1  ZW  *AMC  -  U1*AMW*  ZFACT  -  ZTFETA  * 

2  BARMWCT  )  +  B  2  *  l  AM  G  +ZW  *  YIB  +  BARMWDT  *  ZFACT  )  — B 1  *(  AKET AD T • C TF A CT 

3  * ( EPSF AC T  *  ZFACT  -  ZW/U1J+  ZETADOT  *  C TF AC T/U 1 *( AMW  -  AMC  * 

4  EPSFACT  )+  ZETA  MBARMWDT/Ul  +  EPSFACT  *  Y i 8  )  +  AMETA  ) 

E  COEFFICIENT  CAUCUUATION 

0A(  5 )  =  B2PRIME  *  CTFACT  /Ul  *Ul  *  AMW*  ZTHETA  +  E2*(  Z’TI-'E  T  A*  E  ARM  WO 
IT  -  ZW  *  AMC  +  Ul  *  AMW  *  ZFACT  )  +  Bl*(  ZETA  MEPSFACT  *  AMC 

2  AMW)  +  AMETA  *(  ZW  -U1*EPSFACT  *  ZFACT  )-  AMETADT  *  CTFACT  * 

3  EPSFACT  *  ZTHETA  ) 

F  COEFFICIENT  CAUCUUATICN 

0AI6)  =  Ul*  ZTHETA*  (B2*  AMW  -  AMETA  * 

CO  ICC  J=  1  » 1 29 

100  B(J)  *  C.O 

WRITE  CUTPUT  TAPE  5,3,CTAIU,XKE,UE,B1,E2 
3C  FORMA  T  (///  1 3H  TAIL  CHORD  =  F7.3,1X, 


B 1 *EPSF ACT  ) 


8HK 


iSUB 


SUB  E*  F 9 .5 , IX, 9 KM U 

1 E-  F9 .5* 1 X, 8HB  SUE  1=  F9.5,1X,  EHB  SUE  2=  F9.S  ) 

WRITE  OUTPUT  TAPE  5 , 4 ,82  PR IME » ZE T A, ZE TADOT , AMET A ,  AMETADT 
4CFCRMATI/14H  B  SUB 2  PRIME*  F9.5,1X,  1 1H  Z  XUE  ET A=  F9.5,lX, 

1  15H  Z  SUE  ETA  COT*  F9.S,lX,llH  M  SYB  ETA*  F9.5»1X,16H  M  ISLE  ETA 

2  DOT*  F9.E  ///) 

Cl  *  UFACT  *  XKFACT  «(  Ul  *  AMW  *  (ZETA  -ZTHETA)-  Ull*AN£TA*  ZW) 

C2  =  -B2PRIME  /Ul  *  CTFACT*!  ZW  *AMQ-U1*AMW  *ZFiAOT  -ZTHETA.*EARMWCT) 
C3  =  E2* ( AMC  +ZW* Y  IB  +  BARMWCT*  ZFACT) 

0D4  =  -Bl* (AMET ADT*CT FACT*! EPSF ACT* ZFACT-ZW/U I )+ZETACCT/Ul *CTFACT*< 
1  AMW  -  AMC*EPSFACT)  +  ZET  A* (  EARMWDT /U 1  +EPSF ACT*YI B )  +AMET  A ) 

C  =  D1  +  D2  +C3  +  D  4 
A  (  4  )  *  C 

60F0RMA T  (//  4H  Dl  =  F15.E,  4H  D2=  F15.8,  4H  D3*  F15.8,4H  D4«  F15.8 
1  »  9H  D  T  CT  At.  =  FI  E.8 ,  //  ) 

El  =E2PRIV'  CTFACT  *(U1*AMW* ZTHETA) 

E2  =  B2  *i.-  -  BARMWDT  -ZW*AMQ  +U1  *AMW*  ZFACT  ) 

0E3  *B  1  *  (Zc  in*!  EPSFACT* AM Q  -AMW)+AMETA*(Z W-l!  1*EPSFACT*ZFACT)- 
1  AMETACT*CTFACT»  EPSFACT  *  ZTHETA) 

E  =  El  +  E2  +E2 

70FORMAT  {//  4H  E 1=  F15.8,4H  E2*  F15.8,4H  E3*  F16.8,  9 Hi  E  TOTAL-* 

1  F15.E,  //) 

0C1  *U  F  ACT  *  XKFACT  *<ZW  *AMQ+AMW*  (  ZETACCT*CTFACT-U1  *ZF  ACT  HEARMWCT 
1  «( ZETA-ZTHETA) +  U 1  * AME TA-AMET ADT  *CTFACT *  ZW  ) 

C2  *B2PRIME/U1  *CTF  ACT  » ( AMQ  +ZW*YI8  +ZFACT*  EARMWDT) 

CC3  *  -B2  *YIB  -Bl  *(AM.ETADT/U1  *CTFACT+  ZETACCT/U 1*CTF ACT*  l  YI  E  * 


I 
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1  EPS  FACT  ♦  BARMWCT/U1)  ) 

C  =  Cl  4  C2  +  C3 

WRITE  Ol'T PUT  TAPE  St  8.  Cl.  C2,  C3,  C 
WRITE  OUTPUT  TAPE  5,6,01,02,03,04,0 
WRITE  OUTPUT  TAPE  5,7,  E1,E2,E3,E 

8CF0RMA  T  (///  4H  Cl*  F15.8,  4H  C2*  F15.8,  4H  C3  =  F15.8,9H  C  TOTAL- 
1  F15.8,  //  ) 

RETURN 

END 

SUBROUTINE  POLYMUL  (CRR, CR I , K ,K APPA ,TC AR AT ) 

IMAX  *  50 
DEL  -  C.l 
NUN  *  5 
RATIO  =5.0 
ALTER  *1. C0C001 
EP1  -.CCCOCCOCOOOCCCCCOOCOI 
EP2  -.CCC00CC1 
EP3  *  .OCCCCC  1 
EP4  =.CCCCCCC1 
SR (  1  ) *  -.5 

sim*  c.o 

SR ( 2 ) *  C .5 
SI (2)=  C.O 
SR ( 3 ) *  C.O 
S I  l  3  1  *  0.0 
MODE  *  1 
NNN  =  N+l 

749  IF  (CPP(ll)  750,751,750 

751  CO  753  L*1 , NNN 
M  =  L  +1 
CRR(L)  *  CRR(M) 

753  CONTINUE 

752  N  =  N- 1 
GO  TO  749 

750  CONTINUE 

65  FORMAT  ( 1 C A 6 ) 

DO  5  1=1, 1C 


K  0  N  *  £  K 

I F ( K A PP A ( I  1-KCNJ  2,5,2 

5  CONTINUE 
GO  TC  999 

2  PRINT  65, KAPPA 

66  FORMAT  (2X,2HN'=I4 , 1 2X  ,  5H I  MAX*  1 3 , 1  OX  ,4HNUN=  1 3  ,  1 1  X  ,4HDEL=E8.2  i 

661  6X,6FRAT  IC*E8.2, 4X , 6HAL TER  =  F 1 0. 8/ 

662  2X,4FEP1=E8.2,6X,4HEP2=E8.2,6X,4HEP3=E8.2,6X,4HEP4«E8.2/> 
M0CE=MCDE+1 

67  FORMAT ( 6H  SRI*  E12.5.6H  S I  1=  E12.5,6H  SR2*  E12.5,6H  SI2* 

671  E 1 2 • 5 , 6H  SR3*  E12.5,6H  SI 3*  E12.5) 

NP  1  =  N+  1 

7 1 OFORMA T  (46HCTHE  COEFFICIENTS  OF  THE  GIVEN  PCLYNCMI AL  ARE  /// 

1  1 1H  REAL  PART  ) 

501  WRITE  OUTPUT  TAPE  5,71 

68 10F0RMAT (11 CH  A  B  C  E  £ 

1  F  G  H  /) 

687  CONTINUE 

WRITE  OUTPUT  TAPE  5,  680 , (CRR ( J ) , J*1 , NP 1 ) 

680  FORMAT  ( 1H  8F15. 7  /  ) 

WRITE  OUTPUT  TAPE  5,  682 
682  FORMAT  (16H  IMAGINARY  PART  ) 

WRITE  OUTPUT  TAPE  5*  680 ♦ (CR I ( J ) , J= 1 , NP 1 ) 

502  DO  6  I*  1 , NP 1 
CPR (  I  )  =CRR ( I ) 

6  CPI(I)=CRim 

CALL  C0MAG(CRR(1 ),CRI(11,C1,KE) 

DO  3C5  K* 1 ,N 
GO  TC(9, 609) , MODE 

53  FORMATJ 1H02X4HRCOT5X7HITERANT6X2( 10HREAL  P ART  1 CXl CF INAG  PART10X) 
1,1  OHRCOT  IS  I  N/2X6HN'UMBER4X6HNUMBER7X2<  SH  CF  -ROOT  1  1 X  )  ,2  <  SH  CF  R< 
2z  )  1 1X  )  ,  1 1HA  RADIUS  CF) 

9  PRINT  53 
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6C9 

11 

12 


14 

13 

15 

16 


700 

701 
7C2 

710 

711 

713 

714 

715 
55 

716 
718 


717 

17 

18 

19 


101 

1C2 

103 

104 

105 

106 
107 
1C8 


110 


111 
1  12 
113 


1  =  1 

NPl=N  +  2-K 
NN=NF 1-  1 

IFIK  +  1-N) 13,12,1  1 

CALL  DIVD(-CRR<2),-CRI<2),CRR<1),CRI(  l),XR»mVXI<  1),KE) 
K  1  =  1 
K2=  1 

GO  TC  16 
AR-CPR  (  1  ) 

AI=CRI  (  1  ) 

BR=CPR ( 2  ) 

C I =C  P I ( 3  ) 

K  1  =  1 
K2=1 
M4  =  2 

GOTO  123 
X  R ( 1 )  =  C  E  ARR 
X  I  {  1  )  =  CEAR I 
GO  TC  16 
DO  15  J  =  1 ,3 
XR  t  J )  =  SR ( J  ) 

X I ( J )3SI ( J ) 

K  1  =  1 


M3=  1 
M4=  1 

DO  717  L  =K  1 »  K2 
ZR=XP(L) 

ZI =X  I  (  L  ) 

CALL  PCLYNCM (NP1 , ZR , Z I ,C RR ,C R I ,RR , R I , KE  ) 

FXR(L)=RR 

FXI1 L)=P  I 

CALL  CCPAGIRR.RI .PFAG.KE  ) 

RAC=ALTER*(PPAG/C 1  >*«< 1 • O/FLCATF ( NN) ) 

GO  TC  <715,19, 179), P3 
GO  TC(716,718),MCCE 

FORMAT (  2(16,  4X),  4E20.11,  5X,  E1C.4  ) 

PRINT  55,  K, I,ZR,ZI,RR,RI,RAC 
POLYC=PCLYN 
POLYN=PM  AG 
IFIPPAG 1717,300, 717 
1  =  1+1 

IF(K+1-N) 17,300, 300 

GO  TC  I  10,200  »M  1 

VAL=CEL*  PCLYN 

DBARP=RAD 

D8AR  1=0.0 

K 1  =4 

K2=4 

M 1  =  2 

M3  =  1 

ABARR  =  XR  (  1  1-XR  (3  ) 

ABAR I -XI ( 1 )-XI (3  ) 

BBARR=XR (2  1-XR (3  ) 

B6ARI=XI(2)-XI(3) 

AMIBR  =  XR( 1  )  —  X R (2  ) 

AM  I B I =X I ( 1 )  —  X I  ( 2  ) 

CALL  MVJLT(ABARR,AEARI»BEARR,EBARI  ,TA.TE,K£1  ) 
CALL  MLLTITA ,TB, A N I ER , AM  I B I , DENR ,DEN I ,K E2 ) 

CALL  CCNAGtDENR, DENI ,TA,KE) 

CALL  CCNAGIXR(3) ,XI ( -1 *T4,KE  1 
I F(T/-EP1*T4 11 10 , 1  10, 11  1 

CALL  CERIV(NPUXR<3)  ,XI  (3)  ,CRR,CRI,DR,0I,K50) 
CALL  CCMAG( RR.RI  ,TR,KE) 

CALL  COMAGCDR.CI ,TC.KE) 

IF(TR-EP4*TC) 192 ,171,171 
DELAR=FXRt1)-FXR<3) 

DELA  I=FX  I  (  1  1-FXI  ('31 
DEL0R=FXR(2)-FXR<3) 
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114  DELBI=FXI(2)-FXI (3) 

115  CALL  MULT(eEARR,  e 8 A R I , D E L AR , CEL A  I  , TA , TB , KE5 ) 

116  CALL  MULT(ABARR,  A  E  AR  I  ,  C  E LER  ,  C ELB  I ,  TC  .  T D  ,K;E6  ) 

117  CALL  CIVD(TA-TC, T E-T C , OENR, D EN I , AR, AI.KE7) 

118  CALL  MULT( ABARR, AEARI ,TC,TD,T1 ,T2,KEE) 

119  CALL  MULKEBARR,  EBA R  I  , T A ,TE , T3  ,T 4 , KE 9 ) 

120  CALL  CIVD< T1-T3, T2-T4, OENR, DENI, BR,BI,KE10) 
CR=FXR (3  ) 

CI=FXI  (3  ) 

123  CALL  MULT(ER»BI»ER'»EI,Tl,T2»KEll  ) 

124  CALL  MULT(AR,AI,CR,CI»T2*T4»KE12) 
j*l  TA=Tl-4.C«T3 

132  TB»T2-4 . 0«T4 

CALL  CSCRT (  T  A  ,  TB  ,  TC  ,  TO  ) 

147  T1=-ER+TC 

148  T2*-EI+TD 

149  T3=- ER-T  C 
15C  T  4=- E I-TD 


CALL  CCMAG(T1,T2,TA,KE14) 

CALL  C0MAG(T3,T4,TB,KE15) 

I F ( T  A-T  E  )  154,168,166 
T  A  =  T  E 
T  1  *T  3 
T2*T4 

GO  TC  (  157, 159), M4 
IF(TA)  161, 161, 15  E 
CALL  COM  AG  <2 .0*C  R,2.0*CI,T8,KE16) 

I F (TE-RAD+T A  )  159 ,159,1 8C 

CALL  C I VD ( 2 .0 «CR »  2.0 «CI , T1 , T2 , DBARR, 0BARI,KE17) 

GO  TC  (  161,14)  ,M 4 
XR(4)=XR(3)+C8ARR 
XI (4  )=XI ( 3  )  +  DBAR  I 
TR  =  AESF(XR(4  )  ) 

TI*AESF (XI (4  )  ) 

IFtTR)  167,  167, 14C 
IF ( T I  )  167. 167,169 
I F (TR— T  I  )  164 ,  167,163 
IF (TI-EP2*TR)  165 ,167,167 
XI (4  )  =  C.C 

GO  TC ( 5C3, 504  )  ,MCDE 

F0RMAT(40HCITERANT  ALTERED  TO  BE  PURE  REAL  NUMBER.) 

PRINT  56 
GO  TC  167 

IF(TR-EP2*TI  )  166, 167,167 
XR (4  )=0.0 

GO  TC  (5C5,  167)  ,MCCE 

FORMAT (45HCITERANT  ALTERED  TO  BE  PURE  IMAGINARY  NUMBER.) 

PRINT  57 
GO  TC  7CC 

CALL  0  I V0( T 1 , T2, T  A  , C . C , T 1 , T 2 , K20 ) 

CALL  DIV0(CR,CI, TB  , C. 0 , CR , C I , K2 1  ) 

CALL  MUL T (C R , C I, RAC , C. 0 ,CR , C I , K22) 

GO  TC(5C6,507)  ,MCCE 

FORMAT ( 87H  ITERANT  IS  CUTSICE  CIRCLE  WHICH  BOUNDS  A  ROOT.  INTERP 
5810LATE  ITERANT  TO  ECGE  OF  CIRCLE.) 

5C6  PRINT  58 
GO  TC  159 

CALL  COMAGIABARR, ABARI,T1 ,KR1 ) 

ENA (  1  )  STA(ITA)  EN A ( 3  ) 

I F ( T 1-EP3«T4  )  174 ,1  74,172 
CALL  CCMAG(B8ARR,EEARI,T2 ,KR2) 

ENA ( 2  )  STA(ITA)  ENA ( 3 ) 

IF(T2-EP3«T4)  175, 175,173 
CALL  CCMAGIAMIBR , AM  I  8 1 , T 3, KR3) 

ENA ( 1  )  STA(ITA)  ENA ( 2 )  STA(ITB) 

IF(T3-EP3«T4)  174 ,174,11  1 

ENA(l)  STA(Kl)  STA ( K  2 )  ENA (2 )  STACM3 )  SLJ(17e) 

E  N  A  (  2  )  STA(Kl)  STAIK2)  ENA  (  3  )  STA  (M3  ) 

XRj(K  1  )=XR(K1  )*(1  .0+2 .0«EP3  ) 

XIi(Kl)-XI(Kl)*(1  «C+2.0*EP3) 

GO  TC ( 5C 8 , 5C9 )  .MODE 
6C  FORMAT!  11H0ITERANTS  XIl,6H  AND  XI1.41H  ARE  <TOC  CLOSE  TOGETHER  ALT 
1ER  ITERANT  X  I1’,1H.) 


1  53 
1  54 

155 

156 
168 
1  57 

158 

159 

161 

162 


140 

169 

163 

165 

56 
503 
5C4 

164 

166 

57 
5C5 
167 
180 


58 


507 

171 


172 


173 


174 

175 
178 


STA(ITB) 


STA (ITB) 
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5C8  PRINT  6C,ITA,ITB,ITA 
5C9  GO  TC  7C0 
179  PCLYC=PM AG 
GO  TC  19 

2C0  IF ( PC LYN-VAL ) 201  ,2C1,210 
2C1  VAL=CEL«  POLYN 
2C2  LI  M=  I +NUM 
203  H2=2 

CALL  CERIVINP1  ,XR(4)  ,XI  ( <4 )  ,CRR,CRI ,DR,OI ,K60 ) 
CALL  CCMAG (RR,RI  ,TR, KE) 

CALL  COMAGIDR.DI  .TC,  KE) 

IF(TP-EP4*TC) 192 ,21C,21C 

210  DLT«PAT IO+POLYO/ PCLYN 

211  IF( 1 .O-CLT  )22C, 220,212 

212  CALL  MULT(C8ARR, DEARI,0LT»C.C,08ARR,0EARI,K3C) 
215  L  IM=LIN+ 1 

GO  TC  (510,511), MCCE 

72  FORMAT ( 1 20H0PCL YN CM  I A L  HAS  I  NCR  EASEO  IN  MAGNITUDE 
721RRENT  STEP.  THEREFORE  REOUCE  CURRENT  STEP. 


TOO  MUCH  WITH 


510 

511 
220 
221 
222 

684 


685 


PRINT  72 
GO  TC  161 

GO  TC  (221.231 ),M2 
I F ( I-IMAX)250, 25 0,222 
PRINT  62 
CO  684  J=  1  , NP  1 
CRR(J)  *  1C.C  *  CRR(J) 

ICOUNT  *  ICCUNT  ♦  1 
IF  (  ICCLNT  -  3  )  685,685 ,777 
WRITE  OUTPUT  TAPE  5,686 
6860FORMAT  (  55HCTHE  REVISED  COEFFICIENTS 
1  ///) 

K  =  1 
GO  TO  687 

62  FORMAT! 69HCMAXIMLM  NUMBER  OF  ITERATIONS 
621 P ( Z  )  EY  DELTA.) 


CU 

) 


CF  TEH  HGIVEN  POLYNOMIAL  ARE 


R£»CF£B  WITHOUT  REDUCING 


513 

777 

231 

250 


251 

192 

69 


GO  TC  777 
RETURN 

I F (  I-LIM)250,250,3CC 
DO  2  f 1  L  =  1 ,3 
XR (L  >  =  XR (L+l  ) 

XI ( L )  —  X I (L  +  l  ) 

FXR( L )  =  F  XR ( L ♦ 1 ) 

FXI ( L ) * FX  I  ( L+ 1 ) 

GO  TC  1C1 

GO  TC(514, 300) .MODE 

FORMAT! 19H0FIRST  DERIVATIVE  (  E17.9, E2C.9, 5SH) 
1 RANT  IS  SUFFICIENTLY  CLOSE  TC  ROOT.) 


INDICATES  THAT  ITE 


ARE  AS  FGLL 


514  PRINT  69, DR, Cl 
300  DO  3C2  J=2, NN 

CALL  MULT(ZR,ZI,CRR(J-1 ) ,CR I (J- 1 ) ,TR , T I , K4C ) 

CRR(J)=TR+CRR( J) 

CR  I  ( J  )  =  T  I  +  CR I ( J) 

FORMATC59HOTHE  COEFFICIENTS  OF.  THE  REDUCED  POLYNOMIAL 
631CWS  ) 

ROCTP ( K  )  =ZR 
ROOT  I  (K ) =Z I 
GO  TC(3C3,305),MCDE 
PRINT  63 

FORMAT  (  1  H  6E18.9) 

PRINT  68,  (CRR(  J),  J=1,NN) 

PRINT  68,  (CR I ( J  ),J*1,NN) 

CONTINUE 

GO  TC  (5  15,516  ), MCCE 
PRINT  75 
PRINT  8C 

8CCF0RMATI  1HG42X 1 4HT ABLE  OF  RC0TS//3X4HRCCT8X  1CHREAL  FART1CX10HIMA 
1G  PART  10X  6HPERIDC  10X  26HTIME  TO  CAMP  :T0  HALF  AM.FL 
2  /2X6HNUMEER7X2 (9H  CF  R00T11X)  ) 

'  DO  3C6  I-l.N 

CALL  PCLYNCM  (N+l  , ROOTR  t  I  )  .  RCOT  I  (  I  )  ,  CFR  ,  CP  I »RR»R  1,'KiE  ) 

THALF  =  -.6931471 8C6/R00TR( I )«  TCARAT 
PER  IOC  *  6.2831853C72  /AE SF ( RCOT I ( I ) )  •  TCARAT 
3C6  .WRITE  OUTPUT  TAPE  5,  84,  I,  ROOTR  (I)  *  RCCTI.H  ),  PERI  CC,  THALF 


302 

63 


3C3 

68 

304 

3C5 

515 

516 
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5X'»  F1C.4»5H  SECS) 


75 


84  FORMAT  ( I6,4X,F16 .6 , 3X , F 1 6 .6 , 7X ,F 1 G. 4 , 8H  SECS 
83  FORMAT (16,  4X,  4E2C.11) 

PRINT  75 
FORMAT! 1 H 1 ) 

IMAX  =  50 
NUM  *  5 
CEL  -  C.l 
RATIO  =5. C 
ALTER  =1.  C0C001 
EP1  = .CCCCCCOOOOOCCCCCCCCC1 
EP2  =.CCC00001 
EP3  =.0000001 
EP4  -.0CCCCC01 
SR! 1)=  -.5 
SI l  1)*  0.0 
SR  I  2 ) =  0.5 
S 1 1  2  )=  C.C 
SR  l  3 ) =  0.0 
SI l 3)=  0.0 
tfOOE  =  1 
GO  TC  1 
RETURN 
CONTINUE 
END 

SUeRCUTINE  OER I V!  N , ZR , Z I , CR, Cl , OR , 0 1 , KER ) 

DIMENSION  CR  l  129  )  .CI ( 129  ) 

ENAIC)  STA  (OR)  STA(CI)  STA(RR)  $ T* <  R I )  EN  A(  1  )ST  A  ( K  ER  )  . 

DO  2  J=1  ,N 

CALL  MUL  T( ZR,ZI, RR.RI ,TRR,TR  I, Kl  ) 

CALL  MULTIZR.ZI, CR, C I ,TC R , TO  I , K2 ) 

RSC(Kl)  AJP11L+2)  RSC<  K2  I  AJP  11  )  ENAT2)  ST  Ai(  KER  )  S L J  ( 3  )  . 

DR=TCR*RR 
DI=TCI+RI 
RR=TRR+CR( J) 

R  I  =  TP  I+C I ( J ) 

CONTINUE 
END 

SUERCUTINE  PCL YNOM (tf,  ZR, Z I ,CR,C I ,RR ,R I ,KER ) 

DIMENSION  CR(129  ),CI l 129} 

ENA ( 0 )  STA(RR)  STA(RI) 

DO  2  J«1,N 

CALL  MULT! ZR,ZI, RR.RI, TR,TI,K1  ) 

-  1  }  ENA! 2) 


1 

999 


1 


2 

7 


ENA  ( 1  ) 


STA(KER) 


RSC ( K  1  ) 

1  RR=TR+CR ( J ) 

2  RI*TI+CI(J) 

3  CONTINUE 
ENC 


AJP  l 


STAIKER)  S  L J ( 3 ) 


SUBROUTINE 
CCN (SQ2= 
LCA(XR) 

CSCRT  <XP,XI,YR,YI  ) 
1.4142135624  ). 

LDQIXI)  AJP21L+1) 

LAC(XR) 

CJP21L  +  1 ) 

LCCtXI ) 

STA(A  } 

STQ  l  e  1 

CJP1 ( 1  ) 

STC (C) 

SLJ418)  ♦STA(P)SLJ(£). 

1 

AJP112) 

LCA IB) 

SLJ4 (8) 

+F0VISC2) 

STA (P) 

SSK( X I  ) 

SLJ ( L  +  2  ) 

LAC(P) 

♦  STA ( C  ) 

SLJ15) 

• 

2 

TFSIE) 

SLJ (3) 

STC! S) 

FCV  (B) 

STA! T) 

STA(R)SLJI4) 

3 

STA(S) 

LDA(B) 

FCV! A ) 

STA  (T  ) 

LCA! 1. C) 

STA(R)SLJ(4) 

7 

8 

4 

Y-SQPTFIX) 
SL J (* ) 
LCAtT  ) 

STA  (X) 

FMU (T) 

SLJ ( 7 ) 

FAC (1 .0) 

SLJ4 (8  ) 

♦FAC(R) 

FCV! 2;c) 

SL J4 ( 8 ) 

♦  STA  IT) 

LOA(S) 

SLJ4! 8  ) 

♦FMU(T) 

STA(P) 

5 

LCAIXI) 
SSK! XR) 

FCV  (2.0 
SLJ (6) 

FCV! P) 
LOA(Q) 

STA(C) 

LCC(P) 

A JP2 l L ♦  1 ) 

LAC(C) 

SSXtXI ) 

LDQ(P) 

STA! YR ) 

STC ( Y  I  ) 

StUU+3) 

• 

6 

LCA  IP  ) 

LCQ (Q) 

STA! YR) 

STQ(YI) 

• 

END 

SUeRCUTINE 

LCA(XR) 

CQMAG ( XR*  XI 
LDQIXI ) 

»Z  » KER ) 
AJP21L+  1) 

LAC (XR) 

QJP2 (L+l ) 

LCC! XI ) 

STC(T  ) 

♦THS (T) 

L  LS ( 48 ) 

CJP ( 3  ) 

STC(T) 

FCV ( 7 ) 

STA(H) 

Y-SQRTF(H) 

-FMU(T) 

FMU(H)  FAO(l.O) 

♦EXF7 (141B)SLJ(L+2) 

STA  (H) 

ENC  (1  ) 

SL-JJ3) 

• 

ENC ( 2 ) 

3 

STG(KER) 

END 

STA (Z) 

• 

SUBROUTINE  DIV C  ( XR ,X I  ,  YR,  Y I.,ZR, Z I  , KER  ) 
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2 

1 


1 A 
1 


CALL  PRCC(  XR,XItYR'1-YI,B1,82,PR,PI,DRfDI) 
LCA ( B2 )  A JP  1  ( 1)  ENA ( 3 )  SLJ(S) 

ENA  (2  )  SLJ ( 3 ) 

T*0R«DR+CI«CI 

L  C  A  (  E  1  )  -FDV(B2) 

LCA(PR)  FDV(T) 

LCAIPI)  FCV(T) 

ST  A ( KER  ) 

ENC 

SUBRCUTINE  PUL  T  (XR  ,X  I .  YR,  Y I ,  ZR ,  Z I .  KER  ) 
CALL  PRCClXRfXIt YR,  Y  I , B  1  ,  B2 • PR  ,P 1 ,01 ,C2 ) 

lii\ Hi  -fPB J H }  :f«?  1$  IJibUli 

LCA(PI)  -FPU  (El)  +EXF7 (  1 4 1 B ) SLJ  (  1  ) 
STA(KER)  SLJ  (L  +  2  ) 

STA (KER  ) 


+EXF7(141B)SLJ(2)  STA(Bl) 
-FPU(El)  +EXF7  nU!E)SLg(2) 
-FPU(EI)  ♦  E  X  F  7  ( 1 4  1  B  )  S  L  J  I  2 ) 


STAX  El  ) 
STA(ZR) 
STACZI ) 


EN  A (  1  ) 

ENA (2) 

ENC 

SUERCLTINE  PROC(XR,XI , YR , YI , B 1 , 82 • PR , P I , CR  , Ci ) 
CALL  NCRP(XR,XI,81',AR,AI  ) 

CALL  NCRPXYR.YI,  B2«CR«0I ) 

PRxA  P  »DR-A I »C I 
P  I  =  A  I»CR  +  AR*CI 
ENC 

SUBRCUTINE  NCRP ( A 1 «A2  *  B1 • S 1 • S2) 

SL J ( 1 )  +SEV7 (70000B )  ZRC(C) 

LDA11A+1)  LDQ(Al)  CJP2XL+1)  LCC(Al) 

Q JP2 ( L+ 1  )  LQC ( A 2 )  LCL ( 1 A+l ) +THS ( E ) 

♦  AJPKL  +  2)  STA(Sl)  STA(EI)  SLJ(L  +  5) 

LCA(AI)  FOV(el)  STA(SI)  LDA  ( A2 ) 

END 
END 


ST  A ( ZR  ) 
STA(2I)ENAM) 


+  Z  PC  (  40X30 E  )  ZRC  (  C  ) 
STL(E)  LCCXA21 
SLJU+2)  LDA(E) 

♦  ACC (  1A+2 )  STA(EI) 
F DV ( 81  )  +  STA ( S2  ) 


CC3  1 

RUN  1  CG  26  V  224  FPS  F8U  1.29.63  BELL  KUPBERS 


93.4 

375.  C 

2.393 

1 .0 

288.5 

52.6 

.002376 

32. 

21CC0.0 

.4  772 

8.1 

-.437 

C.1S5 

-.2937 

;S  3  2  C  C  .  C 

224. 

.195 

2.759 

9.3845 

.902 

11.78 

14.08 

29.625 

29.! 

F8U  LCNGITLCINAL 

RUN2  STICK 

FIXED 

M«.209,V=234 

CG24’,  V*  1 

•  1 3VSL 

93.4 

375.0 

2.393 

.95 

232. C 

52.8 

.  CC227  6 

32. ‘ 

22C00.0 

.  4  772 

5.  6 

-.437 

.05  11 

-.380 

S6CCC.C 

234. 

.19 

2.59 

9.3845 

.90 

11.78 

14. C8 

29.375 

29.; 

RUN  3 

C.C-.  24, 

V3 1 39  KTS , 

234  FPS 

3/3/63  LEVEL 

FLITE 

93.4 

375. C 

2.393 

1  .0 

288.5 

52.8 

.002376 

32. 

22000.0 

.4772 

8.1 

-.437 

0.195 

-.380 

96CCC.C 

234. 

’.263 

2.920 

9.3845 

.90 

11.78 

14.08 

29.375 

29.! 
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APPENDIX  III 


SYSTEM  MECHANIZ  ATION  FOR  THE  ANALOG  COMPUTER 
Fig.  1  to  this  Appendix  displays  the  analog  schematic  for  the 
equations  of  motion  of  the  airframe .  Fig  „  2  includes  the  schematic 
for  the  System  1  auto -throttle ,  aircraft  engine ,  and  stick  system 
dynamics  simulations.  Fig.  3  shows  the  System  2  auto-throttle 
schematic . 

Table  I  lists  the  potentiometer  settings  lor  ail  systems » 

Table  II  lists  the  scale  factors  used  in  computing  these  settings „ 
Table  III  describes  the  calculations  used  in  arriving  at  these  values  . 
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TABLE  IT 


SCALE  FACTORS  USED  IN  ANALOG  SIMULATION 

Scale  Factor  Value 


=  .  005  radian  s/volt  =  .  286  degrees/voll 

=  005  raclarts/volt  =  .186  degrees/  vo'iz 

=  .  002  radiar.s/volt  =  .  11?  degrees/voit 

'^Do'v 

=  .005  radian  s/volt  -  .286  de  jrees/  o  * 

=  .002V  «.2%ofV  =  278  KtsA-oli 

c  c 

*X<l 

=  .002  radk-ns/Volr.  -  . 3,15  degrees/  cZ 

°<r 

=  200  lb i /Vo It 

U(Tct 

=  40  lbs/ra.diJiJss  -  volt 

Utce 

=  50  Ibs/kt  -  volt 

o<tcj 

=  12}  e -’•*?. ilf  as  -sec  -  volt 

^t>TCl 

=  2/3  radis/tts  «  sec  -  volt 

°0)TC2- 

2 

=  30  D>s/radl«rv3  ™  sec“  -  volt 

C<DTC3 

n 

=  300  Y:  s/r.ac.'l'iAo  -  sec  -  void 

=  2.3  cfrS/  vo.  t 

&tct 

=  100  It  i/.'o’r  for  System  I  ,  £0  Lhs/voh.  for  System  2 

(Xrcpiiot  —  5 0  L. -./vo 


^Fl 

_  <5  T  V  .. . ,  „  Ts  j>, 

—  &  v 

0<F3 

=  ,  01  ra-ii'ins/vok 

c<DF3 

=  .  1  lb/ ’sec  -  volt 

WpM  F 

=  .73  L>.s/voit 

II 

2? 

. 286  ce cirees/volt 

£ 

t> 

zr 

°r 

ii 

.02  g/volt 

(Xw  = 

.01  ft/ft/sec  -  volt 
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Append ry  !  JJ  Tatle 

P OW  f  £  leu  h'OfiS 

Po^eh'/io^e"^1'  Calculation  d  vJu«_  =  wlesv  , 

Sjaeci  oftews*) 


2. 

3 

4. 

5 

(3 


?CCJ_f(C;  =  -.058 

<*t 

<x rex  £  =  .4 

o(tct 

Ok  £  -  -0/43S 

OCfcOft 

Vg-fts-  R~  -.3433 
2?U*I  *  =  .  6243 

0(u0Ct 

^  _  ..  c  c  / 

.~tfmTTf“  6 


/A) 

in 

i  f\ 
i  A) 

I  m 
J  Ai 


7. 

36 

./  A\ 

8. 

R  —  — 
tf-DX 

641 

1  A! 

7. 

.  A»t <W*  _  . 
//*'  /tfD«. 

990 

.J  M 

10. 

^  £  -  . 

18o 

/  M 

If. 

22t  K&  -  - 
4  rtt 

-2/39 

i  M 

IE. 

0(*  R  _ 
c  b  oO^oCt 

-/49 

Jfl 

13. 

£  = 

tt  ^De-^t 

1  M 

Id. 

(Xc  R  ~ 
*-*>  ^06-*Xt 

.0/47 

1  M 

Rf.-  .1/1 
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